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ABSTRACT 

We present photometric and spectroscopic observations of SN 2007if, an overluminous (My = -20.4), red 
{B-V = 0.16 at Z?-band maximum), slow-rising {tnsn = 24 days) type la supernova (SN la) in a very faint 
{Mg = -14.10) host galaxy. A spectrum at 5 days past Z?-band maximum light is a direct match to the super- 
Chandrasekhar-mass candidate SN la 2003fg, showing Si II and C II at 9000 km s~' . A high signal-to-noise 
co-addition of the SN spectral time series reveals no Na I D absorption, suggesting negligible reddening in 
the host galaxy, and the late-time color evolution has the same slope as the Lira relation for normal SNe la. 
The ejecta appear to be well mixed, with no strong maximum in /-band and a diversity of iron-peak lines 
appearing in near-maximum-light spectra. SN 2007if also displays a plateau in the Si II velocity extending 
as late as +10 days, which we interpret as evidence for an overdense shell in the SN ejecta. We calculate the 
bolometric light curve of the SN and use it and the Si II velocity evolution to constrain the mass of the shell and 
the underlying SN ejecta, and demonstrate that SN 2007if is strongly inconsistent with a Chandrasekhar-mass 
scenario. Within the context of a "tamped detonation" model appropriate for double-degenerate mergers, and 
assuming no host extinction, we estimate the total mass of the system to be 2.4 ± 0.2 Mq, with 1 .6 ± 0. 1 Mq of 
^*Ni and with 0.3-0.5 Mq in the form of an envelope of unburned carbon/oxygen. Our modeling demonstrates 
that the kinematics of shell entrainment provide a more efficient mechanism than incomplete nuclear burning 
for producing the low velocities typical of super-Chandrasekhar-mass SNe la. 

Subject headings: supernovae; general; supernovae: individual SNe 2003fg, 2006gz, 2007if, 2009dc); white 
dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) are of vital importance as 
luminosity distance indicators for measuring the exp ansion 
history of the Universe (Riessefa/. 1998; Perlmutte r et al\ 
IT999I) . They have a small dispersion (^ 0.35 mag) in intrinsic 
peak luminosity, which can be further reduced to 0.16-0.18 
mag by applying a well-known correcti on dependent on the 
width, or decay rate, of the light curve (IPhimps et a/.llT999l: 
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iGuvefg/.l 120071; iJha etal\^2m% . Searches for other lumi- 
nosity correlates with which to derive even more accurate 
luminosity distances from SNe la are underway, with some 
methods d elivering core Hubble residual dispersions as low as 
0.12 m ag (iBailev et a/.l2009UWang et a/.l2009tlFolatelli et al\ 
120091) . SNe la are generally understood to result from the 
thermonuclear explosion of at least one carbon/oxygen white 
dwarf. However, the underlying distribution of SN la progeni- 
tor systems and explosions mechanisms, and the relative rates 
of possible different physical subclasses of SN la events, re- 
main poorly constrained, with potential consequences for the 
average luminosity of events. Any systematic effect which 
may influence the luminosity of different SN la subpopula- 
tions at the level of a few percent has become a cause for con- 
cern for next-gene ration experiments, p articularly redshift- 
dependent effects (iKim ef a/.l 120041: [Under 2004). A better 
understanding of the progenitor systems will place these cor- 
rections on a much firmer conceptual footing and place limits 
on hitherto uncontrolled astrophysical systematics in SN la 
luminosity distance measurements. 

Rare SN la events displaying extreme characteristics or ev- 
idence of unusual physics can often point the way toward 
other, less extreme instances of similar physics which may be 
lurking in the otherwise undifferentiated sample of "normal" 
events. For example, one commonly-invoked rationale for the 
uniformity in pre-correction SN la luminosities is that they 
start with the same amount of fuel, and are triggered by the 
same physics: the SN la progenitor explodes when its mass 
nears the Chandrasekhar limit, Mch = 1 .4 Mq , its mass is com- 
pletely unbound and converted mostly to heavier elements. 
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especially °Ni, the decay of which powers the SN la light 
curve . The single-degenerate scenario (SD; IWhelan & IbenI 
[T973h ensures that the white dwarf slowly approaches Mch 
via accretion from a non-degenerate compa nion. In contras t, 
the double-degenerate scenario (DD; Iben & TutukovlljiH, 
in which two white dwarfs in a binary system merge and ex- 
plode, provides a way for SN la progenitors to exceed Mch 
and to give rise to more luminous events. There are also 
some arguments that single-degenerate, differentially rotat- 
ing white dwarfs with mass exceeding Mca significantly can 
exist ( Yoon & Langer 2005), although the inclusion of baro- 
clinic and magnetohydrodynamic instabilities appears to pre- 
clude this (e.g., Piro..2008 ). There may therefore be a popula- 
tion of SNe la with a distribution of masses greater than Mca, 
with different explosion physics that interferes with luminos- 
ity standardization. The relative rate of such events among 
SNe la in general may also depend on redshift, and unless 
they can be identified or their luminosities accurately cali- 
brated, they need not be common to produce significant biases 
in reconstructions of the dark energy equation of state. 

There are at least four documented examples of over- 
luminous SN la explosions with progenitor mass poten- 
tially exceeding Mca- The first known example of the 
class was S N 2003fg (SNLS-03 D3bb; Howell et al. 2006); 
SN 2 006gz (iHicken ef a/.ll2007h. SN 2007 if (. Akerlof ef q/. | 
20071) . and SN 2009d c fenaka et a/.ll2009l: lYamanaka etM 
2009t ISiiverman et a/.ll20T ^ were discovered later as events 
spectroscopically similar to SN 2003fg. The main evidence 
cited for a very massive progenitor in each of these cases 
was the extremely high luminosity of each of these events, 
and by inference unusually large ""^Ni synthesis. SNe 2003fg, 
and later SNe 2006gz and 2009dc, were noted for being 
overluminous {My ^ -20), with unusually wide light curves 
("stretch" s > 1.1) and C II lines, evidence for unburned car- 
bon, in early or near-maximum spectra. In SN 2003fg, nar- 
row, low-velocity (8000-9000 km s"' ) Si II lines near maxi- 
mum light have been interpreted as evidence for a high grav- 
itational binding energy, further supporting the hypothesis of 
a very massive progenitor; low velocities were also found 
in SN 2009dc. Ejecta velocities inferred from spectra of 
SN 2006gz were closer to those of normal SNe la. 

SN 2007if was dis covered by the Texas Supernova Search 
dAkerlof ef a/.ll2007h in unfiltered ROTSE-IIIb images taken 
on 2007 August 16.3 UT. It was found independently as 
SNF20070825-001 by the Neai'by Supernova Factory (SNfac- 
tory; Aldering et al.\ I2002h . in an ima ge taken in a red 
(RG610) filter with the QUEST-II camera (i Baltav et a/.l2007l) 
on the Palomar Observatory Oschin 1.2 m Schmidt telescope 
("Palomar-QUEST") on 2007 August 25.4. No host galaxy 
was visible in the discovery images; nor in any available 
sky survey images, including the Sloan Digital Sky Survey 
(SDSS), POSS and USNO, making the redshift determina- 
tion and interpretation of early-phase spectra uncertain. An 
optical spectrum we obtaine d with the SuperN ova Integral 
Field Spectrograph (SNIFS; iLantz ef g/.l 120041) on the Uni- 
versity of Hawaii 2.2 m telescope on 2007 August 26.5 UT 
revealed a blue continuum not obviously like a type la super- 
nova. Spectroscopy taken at the Hobby-Eberly Telescope on 
2007 August 29 also failed to identify the nature of the event 
jAkerlof et a/.ll2007l) . However, a later spectrum we obtained 
with the Double Spectrograph on the Hale 5-m telescope at 
Palomar, on 2007 September 6.5 UT, identified SN 2007if 
as a SN la, apparently well before peak. A further SNIFS 



spectrum taken on September 10.5 UT turned out to be an 
unambiguous match to a published spectrum of SN 2003fg 
dHowell et a/.ll200 6b. Cross-correlation of the September 10.5 
SNIFS spectrum with the SNLS spectrum of SN 2003fg sug- 
gested a redshift of 0.070 ± 0.005. 

The faintness of the host, coupled with the unusually large 
luminosity of the supernova, presented a challenge for the de- 
tection of line emission from the host as late as a full year 
after explosion. More recently, however, an optical spectrum 
of the host galaxy was obtained on 2009 August 24.5 with 
the Low Resolution Imaging Spectrograph (LRIS) at the Keck 
Telescope on Mauna Kea, showing Ha and O III A3727 at a 
heliocentric redshift of 0.07416 ± 0.00082. This new redshift 
measurement allows more accurate determination of the SN 
luminosity and the ejecta velocity scale, which in turn enables 
a measurement of the total mass in the explosion. 

In Sections 2 and 3 we present our detailed photometric and 
spectroscopic observations of SN 2007if and its host galaxy. 
In Section 4 we present the bolometric light curve of the SN 
and an estimate of the mass of ^^Ni synthesized in the explo- 
sion, which we find to nominally exceed Mch- In Section 5 
we argue that the red color of the SN, its unusually long (24- 
day) rise time, and the existence of a plateau in the inferred 
photospheric velocity are best explained by the existence of 
an overdense shell in the ejecta, probably caused by the en- 
trainment of an unburned carbon-oxygen envelope. We also 
estimate the total mass ejected in the SN, and the fraction of 
that mass residing in the shell and envelope. Since SN la mass 
estimates are often sensitive to the assumed kinetic energy 
of the explosion, we consider in Section 6 the importance of 
shell structure on arguments associated with mass estimates 
in the literature, and ask whether shell structure may be more 
common in super-Chandrasekhar-mass SN la candidates than 
previously believed. We summarize and conclude in Section 
7. 

2. OBSERVATIONS 

2.1. Optical Spectroscopy 

Observations of SN 2007if were obt ained with SNIFS 
(lAldering eTal] 120021; iLantz et al.\ 12004 '). operated by the 
SNfactory; the observing log is shown in Table [T] SNIFS is 
a fully integrated instrument optimized for automated obser- 
vation of point sources on a diffuse background over the full 
optical window at moderate spectral resolution. It consists of 
a high-throughput wide-band pure-lenslet integral field spec- 
trograph (IFS, "a la TIGER"; Bacon et al. 1995, 2000, 2001), 
a multifilter photometric channel to image the field surround- 
ing the IFS for atmospheric transmission monitoring simul- 
taneous with spectroscopy, and an acquisition/guiding chan- 
nel. The IFS possesses a fully filled 6."4 x 6."4 spectroscopic 
field of view (FOV) subdivided into a grid of 15 x 15 spa- 
tial elements (spaxels), a dual-channel spectrograph covering 
3200-5200 A and 5100-10000 A simultaneously, and an in- 
ternal calibration unit (continuum and arc lamps). SNIFS is 
continuously mounted on the south bent Cassegrain port of 
the University of Hawaii 2.2 m telescope (Mauna Kea) and 
is operated remotely. The SNIFS spectra of SN 2007if were 
reduced using our dedicated data reduction procedur e, similar 
to that presented in Section 4 of iBacon et al.\ ( 120011) . A brief 
discussion on the spectrographic pipeline was presented in 
lAldering et al.\ (|2006), here we outline changes to the pipeline 
since that work, but leave a complete discussion of the reduc- 
tion pipeline to subsequent publications focused on the instru- 
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ment itself. 

After standard CCD preprocessing and subtraction of a low- 
amplitude diffuse-light component, the 225 spectra from the 
individual spaxels of each SNIPS exposure were extracted 
from each blue and red spectrograph exposure, and re-packed 
into two (Xjj, A)-datacubes. This highly specific extraction is 
based upon a detailed optical model of the instrument includ- 
ing interspectrum crosstalk corrections. The datacubes were 
then wavelength-calibrated, using arc lamp exposures ac- 
quired immediately after the science exposures, and spectro- 
spatially flat-fielded, using continuum lamp exposures ob- 
tained during the same night. Cosmic rays were detected and 
corrected using a three-dimensional-filtering scheme applied 
to the datacubes. 

SN and standard star spectra were extracted from each 
(x,y, A)-datacube using three-dimensional point-spread junc- 
tion (PSF) fit photometry over a uniform background (Buton| 
12009). The PSF is modeled semi-analytically as a constrained 
sum of a Gaussian (describing the core) and a Moffat function 
(simulating the wings). The correlations between the different 
shape parameters, as well as their wavelength dependencies, 
have been trained on a set of 300 standard star observations 
in various conditions of seeing and telescope focus between 
2004 and 2007 with SNIPS. This resulted in a chromatic PSF 
empirical model depending only on an effective width (mim- 
icking seeing) and a flattening parameter (for small imag- 
ing defocus and guiding errors). The three-dimensional PSF 
fit takes the atmospheric differential refraction into account 
without resampling. 

During photometric nights, the SN spectra were flux cal- 
ibrated using a flux solution and a mean atmospheric ex- 
tinction derived simultaneously from all spectrophotometric 
standard stars observed during the same night (Buton 2009). 
In non-photometric conditions, an effective attenuation mea- 
surement for each exposure was made using the stars observed 
by the SNIFS multifilter photometric channel. Objects in the 
field (spatially subdivided into five regions each monitoring 
a different wavelength range and treated separately ) were de- 
tected using SExtractor dBertin & Arnoutsll99 6), and the 
produced catalogs matched to (manually-inspected) star cata- 
logs created from deep stacks of the same field. An adapted 
version of the Supernova Legacy Survey photometry code, 
poloka, was then used to determine the so-called photomet- 
ric flux ratio between each exposure and a reference exposure 
chosen from a reference night considered to be photometric. 
A convolution kernel is computed using the matched objects 
between each pair of images in order to make both compara- 
ble flux-wise, meaning that when it is applied to the best qual- 
ity image (the reference one, with best seeing) it homogenizes 
it to the same "photometric frame" as the worst image. The 
photometric flux ratio equals the integral of the convolution 
kernel, an estimation of the relative flux extinction in each fil- 
ter band for the observations on non-photometric nights, and 
is used to correct for attenuation by clouds. Seeing and the 
stability of the atmospheric transmission were assessed using 
quantitative analyses of SNIFS guider video frames acquired 
during our exposures, along with deglitched CFHT Skyprobe 
data. 

To supplement the ANDICAM imaging photometry dis- 
cussed in the next section, we synthesized additional rest- 
frame Bessell BVRI photometry from the SNIFS flux- 
calibrated spectra. The synthetic photometry is shown along- 
side the imaging photometry in Table |2] and in Figures |6] and 

m 



TABLE 1 

Observing log for SNIFS spectra of 
SN 2007IF 



MJD" 


Phase'' 


Exp. Time (s) 


Airmass 


54338.5 


-9.0 


1200 


1.004 


54353.5 


4.9 


1900 


1.003 


54355.5 


6.8 


1230 


1.010 


54358.5 


9.6 


900 


1.013 


54363.5 


14.2 


1740 


1.005 


54365.5 


16.1 


1970 


1.003 


54373.4 


23.5 


1800 


1.007 


54375.5 


25.4 


2000 


1.004 


JHJ / O.J 


28 2 


1 son 


1 ons 


54378.5 


28.2 


1800 


1.025 


54385.4 


34.7 


1800 


1.005 


54390.4 


39.3 


1800 


1.003 


54393.4 


42.1 


2000 


1.027 


54395.4 


43.9 


1800 


1.034 


54395.4 


43.9 


1800 


1.009 


54400.4 


48.6 


1800 


1.017 


54400.4 


48.6 


1800 


1.052 


54403.4 


51.4 


1800 


1.017 


54403.4 


51.4 


1800 


1.052 


54415.3 


62.5 


1800 


1.052 


54420.3 


67.2 


2200 


1.003 


" Observer frame JD ■ 


■ 2400000.5. 





In rest-frame days relative to B-band maximum 
light. 

2.2. ANDICAM BVRI photometry 

Follow-up BVRI photometry of SN 2007if using 
ANDICAM on the CTIO 1.3-m was obtained through 
the Small and Moderate Aperture Research Telescope 
System (SMARTS) Consortium, from 2007 September 10 
through 2007 December 20. Each band was exposed for 
approximately 240 s. Further spectroscopic follow-up was 
also obtained with SNIFS extending through 2007 Novem- 
ber 26. Since the nature and type of SN 2007if were not 
immediately apparent, SMARTS observations began only on 
2007 September 10.5, when the light curve was already in 
decline. Final reference images for subtraction of host galaxy 
light were obtained in the 2008 observing season; although 
the host galaxy is very faint, at late times its contribution 
could still be significant. 

The SMARTS photometry was processed u sing an au- 
tomated pipeHne based on IRAF dTodvl 11993 1). SMARTS 
images were bias-subtracted, overscan-subtracted, and 
flat-fielded using ccdproc. Point sources were de- 
tected (daof ind) and their instrumental magnitudes mea- 
sured (phot), then aperture-corrected to a 6" aperture 
(mkapfile). The images were registered to a stan- 
dard detector coordinate system (xyxymatch, geomap, 
gregister). The resampled final reference images were 
combined (imsum) for added depth in each band, rejecting 
detector defects and cosmic rays by discarding the brightest 
and faintest value at each pixel location in the stack. An abso- 
lute calibration (zeropoint, extinction and color terms ) was es- 
tablish ed on photometric nights from observations of lLandol3 
(119921) standards, fitting a zeropoint and extinction coefficient 
for each night separately as well as a color term constant 
across all nights. The calibration was transferred to the field 
stars for each photometric night separately using the zeropoint 
and extinction but ignoring the color terms, producing magni- 
tudes on a "natural" ANDICAM system which agrees with the 
Landolt system for stars with B-V = V -R = R-I = Q. These 
calibrated magnitudes were then averaged over photometric 
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nights to produce final calibrated ANDICAM magnitudes for 
the field stars. The final reference co-add was normalized 
to each image in turn, with the best-seeing image in each 
pair convolved to match the worst (f itps f,ps f mat ch); the 
(undetected) host galaxy light was then removed by subtrac- 
tion. The positions of all detections of the SN with signal- 
to-noise ratio (S/N) greater than 10 were averaged, and this 
final position used to measure the flux in each image, in or- 
der to achieve more accurate photometry for noisy measure- 
ments at late times. Systematic errors were estimated for the 
resampling, flux normalization and PSF convolution steps by 
measuring the dispersion in the residuals (from a light curve 
constant in time) of the field stars. We estimate an end-to-end 
systematic floor of 1 % on our differential photometry, so that 
our observations of the SN are limited only by sky noise. 

To characterize the ANDICAM system throughput as a 
function of wavelength in each band, we use a product in 
each band of the KPNO filter transmission curve, the quan- 
tum efficiency of a Fairchild 447 CCD as measured by the 
manufacturer (including Mg02 quartz window transmission), 
and the reflectivity of two layers of aluminum. We ad- 
just the central wavelength of each ANDICAM passband us- 
i ng SMARTS observati ons of spectrophotometric standards 
dStritzinger et al\ l2005h . such that our calibrated natural- 
system photometry matches synthetic photometry in the ad- 
justed passbands from the standard star spectra. After these 
steps, the SN's ANDICAM-system EVRI light curve was 
measured by comparison to field stars. The observer-frame 
ANDICAM magnitudes were corrected for Galactic extinc- 
tion u sing E{B-V) = 0.079 (Sch legel. Finkbeiner & DavisI 



119981) and Ry = 3.1 (Cardelli ef a/. 1989|), then^ 



corrected (*Nu gent. Kim. & Perlmutted 12003) to rest-frame 
Bessell BVRI bandpasses (Bessell 1990), using the adjusted 
ANDICAM filter bandpasses and the SNIFS spectral time se- 
ries. 

2.3. ANDICAM J photometry 

The SMARTS observations were taken with ANDICAM in 
dual CCD-IR mode, so that a J light curve was obtained si- 
multaneously with the BVRI light curve. Four dithered 50 s 
7-band exposures were taken during each 240 s optical-band 
exposure, for a total J exposure time of 800 s each night. 

The y-band images were reduced using IRAF. Separate bias 
and overscan levels for each of the four amplifiers of the de- 
tector were subtracted from the ANDICAM images, which 
were then flat-fielded using a superflat made from dome flat 
images taken nightly. Bad pixels were identified as 3cr-deviant 
pixels in the superflat. Sources in the field were generally de- 
tected at low S/N, precluding automatic registration of sin- 
gle dithered exposures. We therefore performed an initial 
co-add of the exposures at each dither position without fur- 
ther registration; we estimate a positional error of approxi- 
mately 0."5, resulting in a modest loss in S/N. Sources were 
then detected in the co-add at each dither position, and the 
single-dither co-adds were registered to each other to pro- 
duce a final co-add for the night. Host galaxy subtraction 
was neglected for 7-band, since registration was difficult and 
since the high sky noise in these measurements dominates 
over any systematic error from unsubtracted light from such 
a faint host. The field star magnitudes were calibrated us- 
ing standard star observations on nights reported as photo- 
metric by the SMARTS queue observer, with extinction cor- 
rections based on a seasonal average 7-band extinction coef- 
ficient (O.I mag airmass"': iFrogel..l998.) for CTIO. The SN 



magnitudes were then measured via differential photometry 
using three suitable non-variable stars in the FOV, similar to 
the process used in reducing the optical data. 

While we have no near-infrared spectra of SN 2007if with 
which to perform A'-corrections on the 7-band magnitudes, we 
use th e spectra of the 1991T-Hke SN 1999ee (Hamuvefa]] 
l2002h as a reasonable approximation. The transfer function 
used was the product of the quantum efficiency of the Rock- 
well HgCdTe "Hawaii" array and the 7-band filter transmis- 
sion as measured by the manufacturer 

Despite this somewhat approximate treatment, we expect 
based on residuals for the field stars that our systematic errors 
are of order 10% or less, in general comparable to or less than 
the sky noise. These are included in our stated error bars. 

2.4. Additional optical photometry 

To further constrain the SN's multi-band behavior near 
maximum light, the ANDICAM data were supplemented with 
observations near maximum from the broad-band RG610 
data collected by the Palomar-QUEST supernova search, and 
with additional V-band photometric observations taken with 
the SNIFS imaging channel when adjusting the telescope's 
pointing to place the SN into the FOV of the SNIFS inte- 
gral field unit. Stars in the RG610 images of the SN were 
matched to corresponding /?-band and /-band images taken 
with ANDICAM, and their R magnitudes were corrected to a 
natural Palomar-QUEST RG610 system using a linear color 
correction, RG61Q = R + c(R-I) (so that RG610 = for a star 
with/?-/ = 0). The RG61Q magnitudes of the SN in the search 
observations were established via comparison with several lo- 
cal field stars. Finally, the SN RG610 magnitudes were K- 
corrected to ANDICAM R, using /T-corrections calculated 
from the measured transmission curves for the ANDICAM 
R and Palomar-QUEST RG610 appHed to SNIFS spectra on 
either side of maximum light. The process has an estimated 
error floor of about 0.05 mag. A similar procedure was fol- 
lowed for the SNIFS V-band photometry. The RG610 and 
SNIFS V photometry have been merged, respectively, with 
the ANDICAM R and V lightcurve points in the analysis to 
follow. 

FinaUy, lAkerlof ef a/.l (l2007h and Yuan et al. (2010, 
in preparation) present unfitered photometry obtained with 
ROTSE -IIIb that is valuable in constraining the rise time (see 
Sectionl4!2b. 



2.5. Host galaxy observations 

The host of SN 2007if is not visible in pre-explosion sur- 
vey images, and no host spectroscopic features were detected 
in SNIFS spectra obtained while the SN was being actively 
followed. On 2009 August 24.6 UTC, imaging obt ained with 
the L ow Resolution Imaging Spectrograph (LRJS: lOke et alJ 
1 1 9951) on the Keck-I 10 m telescope revealed an extended 
source coincidence with the original SN location. Five ex- 
posures each of 100 s duration were obtained using the LRIS 
blue-side camera equipped with a g-band filter. The expo- 
sures were dithered to allow rejection of cosmetic defects. 
The raw images were overscan-subtracted, flat-fielded us- 
ing dome fl ats, astrometrically calibrated using WCS Tools 
dMinkl I2OO6 ) and then coadded using SWarp (.Bertin et alj 
I2OO2IK The final ima ge, shown in Figure [l] w as segmented 
using SExtractor dBertin & Arnoutsl[T996 l) and then flux 
calibrated using stars overlapping with SDSS-I, matched us- 
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TABLE 2 

Rl:ST-l'RAMI;BVR/,7 LIGH T ri RVli Ol- SN 2()07ll-- 



MJD" 


Phase'' 


B 


V 


R 


/ 


J 


Instrument 


54337.3 


-10.3 






17.78±0.10 






QUEST 


54338.5 


-9.2 




n. 61 ±0.02 








SNIFS-P 


54338.5 


-9.2 


17.54 ±0.04 


17.62 ±0.06 


17.53 ±0.04 


17.64 ±0.06 




SNIFS-S 


54346.3 


-2.0 






17.15 ±0.07 






QUEST 


54353.4 


4.6 






17.29 ±0.05 






QUEST 


54353.5 


4.7 


17.45 ±0.04 


17.29 ±0.05 


17.27 ±0.04 


17.46 ±0.05 




SNlFS-S 


54354.4 


5.5 


17.49 ±0.02 


17.28 ±0.02 


17.28 ±0.01 


17.40 ±0.02 




SMARTS 


54354.4 


5.6 






17.19±0.11 






QUEST 


54355.5 


6.6 




17.29 ±0.02 








SNIFS-P 


54355.5 


6.6 


17.49 ±0.04 


17.24 ±0.04 


17.21 ±0.03 


17.38 ±0.05 




SNIFS-S 


54356.3 


7.3 


17.57 ±0.03 


17.27 ±0.02 


17.29 ±0.01 


17.43 ±0.02 


18.41 ±0.17 


SMARTS 


54358.2 


9.2 


17.69 ±0.04 


17.33 ±0.03 


17.34 ±0.01 


17.49 ±0.02 


18.43 ±0.17 


SMARTS 


54358.5 


9.4 




17.36 ±0.02 








SNIFS-P 


54358.5 


9.4 


17.66 ±0.09 


17. 19 ±0.09 


17.18±0.05 


17.34±0.10 




SNIFS-S 


54361.3 


12.0 


17.96 ±0.02 


17.50 ±0.03 


17.47 ±0.01 


17.58 ±0.03 




SMARTS 


54363.2 


13.8 


18.17±0.03 


17.55 ±0.03 


17.47 ±0.04 


17.60 ±0.03 


19.16±0.14 


SMARTS 


54363.5 


14.0 




17.59 ±0.04 








SNIFS-P 


54363.5 


14.0 


18.14 ±0.07 


17.55 ±0.06 


17.48 ±0.04 


17.59 ±0.07 




SNIFS-S 


54365.5 


15.9 


18.29 ±0.09 


17.55 ±0.06 


17.46 ±0.04 


17.52 ±0.07 




SNIFS-S 


54366.3 


16.6 


18.52 ±0.06 


17.71 ±0.05 


17.62 ±0.05 


17.68 ±0.11 




SMARTS 


54372.3 


22.2 


19.05 ±0.32 




17.68 ±0.23 


17.55 ±0.29 


19.08±0.13 


SMARTS 


54373.4 


23.3 


19.02±0.12 


18.06 ±0.08 


17.76 ±0.05 


17.63 ±0.06 




SNIFS-S 


54375.2 


24.9 


19.38 ±0.08 


18.22 ±0.04 


17.85 ±0.02 


17.73 ±0.04 


19. 17 ±0.20 


SMARTS 


54375.4 


25.2 


19.15±0.11 


18.16 ±0.06 


17.85 ±0.04 


17.73 ±0.05 




SNIFS-S 


54378.2 


27.8 


19.49 ±0.06 


18.37 ±0.03 


17.98 ±0.02 


17.85 ±0.03 




SMARTS 


54378.5 


28.0 


19.38±0.17 


18.33 ±0.09 


17.99 ±0.04 


17.83 ±0.07 




SNIFS-S 


54382.3 


31.5 


19.65 ±0.07 


18.46 ±0.03 


18.11 ±0.02 


17.85 ±0.03 




SMARTS 


54385.4 


34.4 


19.57±0.14 


18.48 ±0.07 


18.11 ±0.04 


17.87 ±0.06 




SNIFS-S 


54388.1 


37.0 


19.81 ±0.07 


18.62 ±0.05 


18.25 ±0.03 


17.99 ±0.03 


19.16±0.10 


SMARTS 


54388.2 


37.1 




18.59 ±0.07 








SNlFS-P 


54388.4 


37.3 




18.64 ±0.07 








SNlFS-P 


54388.4 


37.3 


19.84 ±0.43 


18.69 ±0.22 


18.26±0.10 


17.98 ±0.14 




SNIFS-S 


54390.4 


39.1 




18.71 ±0.05 








SNlFS-P 


54390.4 


39.1 


19.70 ±0.1 5 


18.66 ±0.09 


18.31 ±0.05 


18.05 ±0.07 




SNIFS-S 


54392.2 


40.8 


19.88±0.10 


18.81 ±0.04 


18.38 ±0.03 


18.14±0.04 




SMARTS 


54393.3 


41.8 




18.81 ±0.05 








SNIFS-P 


54393.4 


41.9 


19.79 ±0.26 


18.77 ±0.14 


18.39 ±0.07 


18.13±0.10 




SNIFS-S 


54395.1 


43.5 


19.88±0.23 






18.45 ±0.13 




SMARTS 


54395.3 


43.7 


19.83 ±0.18 


18.82±0.10 


18.48 ±0.06 


18.22 ±0.07 




SNIFS-S 


54400.2 


48.2 


20.49 ±0.49 




18.68 ±0.25 


18.38±0.25 


20.55 ±0.27 


SMARTS 


54400.4 


48.4 




19.06 ±0.09 








SNIFS-P 


54400.4 


48.4 


19.74 ±0.51 


19.19±0.21 


18.90±0.11 


18.71 ±0.13 




SNlFS-S 


54403.4 


51.2 




19.09 ±0.04 








SNIFS-P 


54403.4 


51.2 


20.07 ±0.51 


19.14±0.22 


18.83±0.13 


18.60±0.16 




SNIFS-S 


54404.2 


52.0 


20.04 ±0.10 


19.18±0.06 


18.75 ±0.04 


18.61 ±0.06 


20.43 ± 0.26 


SMARTS 


54408.2 


55.6 


20.20 ±0.10 


19.26 ±0.05 


18.94 ±0.04 


18.61 ±0.05 


20.77 ± 0.33 


SMARTS 


54412.3 


59.5 


20.20 ±0.10 


19.36 ±0.05 


19.08 ±0.04 


18.97 ±0.07 


20.99 ± 0.45 


SMARTS 


54415.3 


62.3 




19.28 ±0.05 








SNIFS-P 


54415.3 


62.3 


20.31 ±0.33 


19.43 ±0.20 


19.19±0.12 


18.98±0.17 




SNIFS-S 


54418.1 


64.9 


20.28 ±0.10 


19.50 ±0.05 


19.23 ±0.04 


19.22 ±0.07 




SMARTS 


54420.3 


66.9 




19.50 ±0.04 








SNIFS-P 


54420.3 


67.0 


20.25 ±0.22 


19.42 ±0.14 


19.25 ±0.09 


19.09 ±0.12 




SNIFS-S 


54422.1 


68.6 


20.16±0.14 


19.50 ±0.08 


19.41 ±0.06 


19.26 ±0.09 




SMARTS 


54430.4 


76.4 




19.49 ±0.10 








SNIFS-P 


54431.1 


77.0 


20.40 ±0.27 


19.60 ±0.14 


19.65 ±0.14 


19.38±0.13 




SMARTS 


54435.1 


80.7 


20.67 ±0.17 


19.76 ±0.07 


19.69 ±0.07 


19.63±0.15 




SMARTS 


54438.1 


83.5 


20.70±0.15 


19.73 ±0.06 


19.80 ±0.06 


19.74±0.12 




SMARTS 


54442.1 


87.2 


20.53 ±0.14 


19.82 ±0.07 


19.96 ±0.08 


20.23 ±0.23 




SMARTS 


54445.1 


90.0 


20.84 ±0.18 


20.00 ±0.08 


19.81 ±0.07 


19.81 ±0.14 




SMARTS 


54448.1 


92.8 


20.87 ±0.1 8 


19.92 ±0.07 


19.98 ±0.07 


19.88±0.14 




SMARTS 


54451.1 


95.6 


20.84 ±0.33 


19.82±0.12 


20.00 ±0.14 


19.79 ±0.18 




SMARTS 



" Observer frame ID -2400000.5. 
In rest-frame days relative to i5-band maximum light. 
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Fig. 1 . — Image of the host galaxy of SN 2007if as co-added from 
Keck LRIS imaging (north is up, east is left). The FOV is 43" x30". 

ing WE S I >0. The resulting host location is a = 1 : 1 0:5 1 .4 1 2, 
S = 15:27:39.57 and the magnitude based on SExtractor 
mag_iso_cor after correcting for Galactic extinction of 
Ag = 0.34 is ^ = 22.89 ± 0.04. The host major and minor axes 
are 1 '.'5 and 1 "2, respectively, uncorrected for seeing of 0."76 
FWHM. 

Upon detection of the host, LRIS was reconfigured for 
spectroscopic observations intended to measure the host red- 
shift. On the blue side the 600 1 mm"' grism blazed at 4000 A 
was employed, covering 3500-5600 A, while on the red side 
the 900 1 mm"' grating blazed at 5500 A was used to cover 
5440-7640 A. The D5600 dichroic was used to direct light 
to the appropriate channel. A slit width of 1" was used, giv- 
ing resolutions of A/AA ~ 1200 and ^ 1650 for the blue and 
red sides, respectively. Four exposures of 900 s each were 
obtained. These were overscan-subtracted, flat-fielded using 
internal quar tz flats, cosmic-ray rejected using LAcospec 
dvan DokkumI 120011) . and then wavelength-calibrated using 
arc and night sky lines. The reduced spectrum revealed weak 
emission lines from Ha and O II A3727. The resulting he- 
Uocentric redshift is z = 0.07416 ± 0.00082, and the result- 
ing host luminosity is Mg = -14. 10 ± 0.07 after correction for 
Galactic extinction. 

3. DISCUSSION 

Our extensive photometric and spectroscopic dataset allows 
us to examine several unique features of SN 2007if. We 
begin with comparison to the spectrum of the prototype — 
SN 2003fg — as well as of SN 2006gz, focusing in particular 
on the evidence for unburned carbon in such systems. The ve- 
locity evolution is examined, especially in light of the low ve- 
locities previously found for SN 2003fg and SN 2009dc. We 
then examine the photometric properties and the influence of 
dust extinction, in preparation for determining the mass of the 
SN 2007if progenitor. 

3.1. Spectral features and velocity evolution 

Figure |2] displays selected SNIFS spectra of SN 2007if, 
along with selected published spectra from SN 2006gz and 
the only spectrum of SN 2003fg available. The SNIFS spec- 
tra of SN 2007if are supplemented by the +1 d spectrum ob- 
tained at the Palomar Hale 5.1-m telescope with the Double 

'^ |http : //nvogre .astro . Washington . edu : 8080/wesix/| 



Spectrogr aph. The spectral t ime-series of SN 2006gz pub- 
lished by iHicken et al] (l2007h extends from -14 d to +1 1 d 
with respect to maximum — the earliest spectrum from that 
event had more pronounced absorption features clearly iden- 
tifying it as a SN la, and also a prominent C II A6580 ab- 
sorption initially attributed to a low velocity Si II component 
jPrieto & Depov 2006). In contrast, the earliest SNIFS spec- 
trum of SN 2007if is relatively featureless, and without prior 
knowledge of the host redshift, not obviously from a SN la. In 
general, the absorption features common to both SN 2006gz 
and SN 2007if appear to be broader, stronger, and more highly 
blue-shifted in the former than in the latter — note the faster 
Si II absorption around 6100-6200 A in SN 2006gz. Around 
maximum, the spectrum of SN 2007if seems to more closely 
resemble that of SN 2003fg than it does SN 2006gz, in partic- 
ular at the double-notch absorptions around 4000^200 A, but 
also note the double absorption at 3700-3900 A missing from 
SN 2007if. Comparison at much later phases to SN 2006gz is 
not possible since that time-series ends much earlier than that 
from SNIFS. The detailed evolution at these phases is rather 
slow, but these spectra appear quite similar to those of a nor- 
mal SN la except for the lack of emission around 4000 A. 

In Figure [3] we present the detailed evolution of the SNIFS 
spectra in the region of the Si II A6355 feature. In the ear- 
liest spectrum, the feature with the blue absorption edge at 
6360 A may be a sign of C II A6580 at a velocity of ^ 
10,000 km s"' — suggesting the presence of unburned car- 
bon as in SNe 2006gz and 2009dc, though the line is much 
weaker. Due to the redshift of SN 2007if, this region overlaps 
telluric B-band absorption, but we have confirmed the accu- 
racy of our telluric correction and believe the notch to be a 
real SN feature. The evidence for C II at this wavelength in 
SN 2003fg is more circumstantial owing to that spectrum's 
lower S/N at these wavelengths and the putative identifica- 
tion of C II A4267 (but see below). Immediately blueward is 
the expected position of Si II A6355, at which an absorption 
"slump" is detectable. The sloping continuum and proximity 
to C II make the line profile difficult to measure directly; we 
factor out the continuum in our measurements of the absorp- 
tion minimum described below, and include related uncertain- 
ties in our error bars. In any case, the shallow inflection near 
9,000 km s"' and the blue edge near 1 1,000 km s"' (both rela- 
tive to 6355 A), coincident with the edge of the Si II feature in 
the post-maximum spectra, suggest that the photosphere may 
have receded to the post-maximum velocity at a phase of -9 d. 
The relative robustness of this edge across phase may suggest 
that the layer of Si-peak elements production extends only to 
11,000 km s"'. 

Figure]?] compares the evolution of the position of the Si II 
A6355 absorption minimum in SNe 2007if, 2006gz, 2009dc, 
and 2003fg. At least in the case of SN 2007if, by the time of 
the next SNIFS observation after day +16 (on day +24), the 
SN spectrum has evolved too much to associate Si II A6355 
with any given spectral feature without peril. We exclude the 
DBSP spectrum from this analysis, since it was obtained at 
lower resolution and in poor conditions; however a measure- 
ment from the profile is consistent with the SNIFS observa- 
tions a few days later. The measurements of Si II A6355 on 
SNIFS spectra were made as follows: bins in each spectrum 
immediately to the right and left of the line feature were used 
to fit a linear local continuum, Fx^cont = a+bX. The spectrum in 
the region of the line was then fitted using F\ = ^A^cont x (/"(A) 
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9 d (SN2006gz) 



2 d (SN2006gz) 
+1 d (DBSP) 
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Fig . 2. — Selected rest-frame SNIPS spectra of SN 2007if (black curves), int erleaved with obser vations of SNe 2006gz (red. iHicken et al\ 
120071) and 2003 fg (blue, the prototypical super-Chandrasekhar SN la candidate. iHowell et a/.ll200H) . A spectrum obtained just after maximum 
with the Double Spectrograph (DBSP, violet) on the 5.1 m Palomar Hale telescope is also included. Spectra of SN 2006gz are selected based 
on proximity in rest-frame phase for comparison. The -9 d spectrum of SN 2006gz has more developed absorption features than SN 2007if 
at the same phase, making SN 2007if more difficult to recognize as a SN la. SN 2007if is in some respects a better match to SN 2003fg than 
SN 2006gz, but not in every case. Note especially the more blueshifted absorptions of SN 2006gz (especially Si II around 6100 A). All spectra 
shown have been binned to 1000 km s"' and later phase SNIPS spectra combined for presentation purposes. 
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Fig. 3. — Detailed evolution of SN 2007if spectra around the Si II 
A6355 feature. 
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Fig. 4. — Comparison of tlie evolution of the Si II absorption feature 
in SN 2007if (black circles) with SNe 2006gz (red squares), 2009dc 
(green diamonds) and 2003fg (orange triangle). The mean evolu- 
tion o f Si II velocity in normal "LVG" subclass SNe la jBenetti et al.\ 
120051) is presented as the dashed line. 
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FIG. 5.— SYNOW fit to the SNIPS day +5 spectrum. 



where the absorption line was modeled as a skewed Gaussian, 



(j)(X)=l-c[l+diX-Xor]e 



i-\„-(X-\of/2cr- 



(1) 



The line velocity was determined by solving d(f)/dX = via 
the Newton-Raphson method. The error bars on the proce- 
dure were determined through a bootstrap Monte Carlo in 
two stages. In the first stage, values of a and b representing 
possible continua ^A.cont were sampled using the covariance 
matrix of the local continuum fit; for each candidate contin- 
uum, the line profile (piX) was then refit holding a and b fixed, 
and values of c, d, Xq and a were sampled using the covari- 
ance matrix of the fit to (f>(X). The final velocity values and 
their eiTors were determined as the mean and standard de- 
viation of the distribution of absorption minimum velocities 
thus generated. We can see in Figure |4] the measurements of 
the velocity and their error with this method for each spec- 
trum of the study. Clearly, the average measured characteris- 
tic ejecta velocities of SNe 2007if and 2009dc are much more 
consistent with that measured for SN 2003fg, the prototypical 
super-Chandrasekhar event. The shallow slope of the velocity 
evolution (v = 34 ± 15 km s ~' day"') is charact eristic of the 
low-velocity-gradient (LVG: [Benetti et a]1l2005h subclass. 

To further characterize the SN ejecta, we have analyzed 
the -1-5 day SNIFS spectrum using the automate d SN spec- 
troscopic direct analysis program, synapps (Thoma s et al.\ 
2009). This code incorp orates the familiar SYNOW- style 
(Jeffery & Branch 1990; Branch. Baron & Jefferyll200l pa- 
rameterized model as part of an objective function in a mul- 
tidimensional nonlinear optimization problem. From a good 
fit the presence of given ions may be ruled out, and for ions 
that are positively identified, corresponding intervals of ejecta 
velocity are constrained. Such a parameterized approach is in- 
dispensable in the analysis of such particularly unusual events 
as SN 2007if for which no reasonable detailed theoretical pre- 
diction otherwise exists, and will continue to play a vital role 
as long as new types of transients are discovered. 

A fit to the day -1-5 spectrum appears in Figure|5] The entire 
SNIFS wavelength range is included in the fit. Familiar SN la 
ion species Mg II, Si II, S II, Ca II, Fe II and Fe III produce un- 
ambiguous spectroscopic absorption signatures. Blanketing 
by iron-peak elements is invoked to suppress UV flux, but the 
optimized model may overcompensate for the tendency of the 
Sobolev approximation to underblanket here. The Si II and 
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S II features are rather weak, a trait in common with the so- 
called overluminous or SN 1991T-like SNe la. The Ca II H&K 
absorption at 3800 A is extremely sharp and weak compared 
to normal events, as was noted in the case of SN 2003fg. The 
highly blended regions centered on 4500 A and 5000 A are ac- 
tually dominated by Fe III, again harkening to SN 1991T-like 
events. The relatively weak IME signatures and the preva- 
lence of Fe III indicate high temperatures in the ejecta; in 
SN 199 IT this has been attributed to the presence of ^^Ni in 
the outer layers of the ejecta (.Mazzali et a/.iri995l) . We also 
detect a weak signature of C II at around 8000 km s"', i.e., at 
velocities comparable to the Si II feature. 

Of particular interest is the pair of sharp notches at 4030 A 
and 4130 A. In iHowell et al\ ( 12006) . the bluer notch is at- 
tributed to Si II, and the redder to C II. In the fit presented 
in that work, a strong C II A6580 absorption was predicted, 
but the mismatch between the prediction and observations 
was attributed mainly to the poor S/N in those spectra. The 
S/N in the SNIFS specti'um is much higher, but the C II 
A6580 is very weak . While NLTE effects may be in play 
jThomas et al\ l2007h . we considered further alternative ex- 
planations for the 4130 A feature. An interesting possibility is 
Cr II absorption — the inclusion of this ion also helps produce 
the strong emission at 4600 A, and contributes to the absorp- 
tion at 3280 A, along with other iron-peak elements. Cr II is 
not without precedent in SN la spectra, though its appearance 
in more no rmal events seems limited to postmaximum epochs 
dBranch et a/.ll2008l) . However, detailed reconstruction of the 
4000-4130 A region remains difficult — in particular, Ni II is 
invoked to explain some of the UV line blanketing, but pro- 
duces an extra notch at 3950 A. 
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Fig. 6.— Rest-frame Bessell BVRIJ light curve for SN 2007if from 
ANDICAM+SNIFS-^QUEST. Upright triangles: QUEST RG6\Q, 
A'-corrected to rest-frame Bessell R. Diamonds: SNIFS V from 
photometric imaging channel, A^-corrected to rest-frame Bessell V. 
Inverted triangles: Rest-frame Bessell BVRI magnitudes synthe- 
sized from SNIFS flux-calibrated spectroscopy. Circles: ANDICAM 
BVRIJ, A^-corrected to the respective rest-fra me Be ssell filters. 
Squares: ROTSE-III points from Akerlof eta!} ( 120071) . A fourth- 
order polynomial fit to the y-band data is also shown. 



3.2. Maximum-light behavior, colors and extinction 

The rest-frame Bessell BVRIJ light curve of SN 2007if is 
given in Table |2] and in Figure |6] The color evolution is 
shown in Figure |7] An accurate estimate of the luminosity 
and ^^Ni mass for SN 2007if requires some care, due in part 
to the sparse multi-band light curve coverage around maxi- 
mum light. In this section we explore two different methods 
of interpolating the light curves to maximum, derive the ob- 
served color of the SN, and comment briefly on reddening by 
dust in the host galaxy. 

Fitting the multi-band light curve of SN 2007if with exist- 
ing SN la templates trained on normal SNe la is of course 
problematic. As an example, the SALT2 model (iGuv et all 
I20b7h fit to all four bands is poor (xl = 930.3/122 = 7.6), 
above the threshold used by the SNL S to photometr ically 
screen events as "probable" SNe la (Sul livan et al.\206^ . The 
event would therefore have been deprioritized for spectro- 
scopic follow-up by SNLS and might not have been studied 
further The fit is especially poor in R and /, which lack the 
distinct second maximum typical of nor mal SNe la; such be- 
havior also was note d in SN 2006gz bvlHicken et al.\ (l2007l) 
and in SN 2009dc bv lYamanaka et al.\ (120091) . 

We can nevertheless attempt to use SALT2 simply as a 
method for interpolating the shape of the light curve around 
maximum light, taking the fitted parameters xi and 

c provisionally. To minimize the impact of details of the 
SALT2 spectral model on the outcome, we use SALT2 in 
the rest frame, include SALT2 light curve model errors 
in the fitting, and we fit BVR bands only, excluding the 
clearly discrepant /. This procedure gives a more reason- 
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Fig. 7. — Rest-frame Bessell color evolution for SN 2007if from 
ANDICAM (circles) and SNIFS spectrophotometry (inverted trian- 
gles). A fit to the Lira relation with £(S-V) = 0.18is shown. 
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able /v= 144.8/93 = 1 .6, with a date of B-band maximum 
around MJD 54348.4 (2007 September 6), a stretch of 1.15 
{x\ = 1.83, consistent with SN 2003fg) and a very red color 
(c = 0.24). The best-fit magnitudes at B-band maximum are 
Mb = 17.34 ± 0.04, my = 17.18 ± 0.05, mr = 17.14 ± 0.06 
(statistical + /T-correction errors added in quadrature), giv- 
ing B - V = 0.16. The value of niR so derived is also con- 
sistent with the observed unfiltered magnitude at maximum 
(calibrated to USNO-B 1.0 R Monet, D. G. et al. 2003) ob- 
served by ROTSE-IIIb ( Akerlof ef a/.il2007rt and the QUEST- 
II RG610 point at that date. 

As a cross-check, we fit a cubic polynomial to the rest- 
frame BVR points before MJD 54374.0 (three weeks after 
the SALT2 date of B-band maximum). This gives maxi- 
mum magnitudes niB = 17.33 ± 0.03, my = 17.25 ±0.03, and 
niR = 17.19 ±0.02, with the date of V-band maximum consis- 
tent with the SALT2 fits to within ±2 days. The uncertainties 
on the date of maximum are larger with this approach, since 
each band can vary independently; however, the results are 
roughly consistent with SALT2. We take the results of the 
SALT2 fit as our fiducial values in further discussion — in 
particular the phase of maximum light at MJD 54348.4, in 
which case the SMARTS observations begin at day +5 after 
maximum light. 

Both of the above extrapolations imply that SN 2007if is 
redder than the typical SN la at maximum light. To as- 
sess the im pact of dust extinction, we use the Lira relation 
(iPhillip s et al. 1999; Folatelli et al. 2009) and t he equivalent 
width of Na I D absorption (Turatto. Benetti & Cappellard 
12002) ("TBC relation"). Two possible slopes for the TEC re- 
lation are given, one shallow (E(B-V) ^ 0.5 x EWCNa I D)) 
and one steep (E(B-V) - 0. 15 x ^^(Na I D)). 

The B-y color evolution of SN 2007if (see Figure E) has 
a slope compatible with the the Lira relation (P hillips et al.\ 
II999) with a fitted color excess of E{B -V\osi = 0.18 ± 
0.04 (stat) ± 0.06 (sys). Assuming the Lira relation holds 
(which for this unusual SN is by no means clear), the intrinsic 
color of the SN is B - V = -0.02. 

A noise-weighted co-add of the SNIPS spectra reveals no 
Na I D absorption at the redshift of the host. To derive an 
upper limit on ^^(Na I D) while allowing for uncertainty in 
the redshift, we construct the probability surface over equiv- 
alent width and redshift for a joint fit of the Na I D line 
profile within the appropriate wavelength range for all ob- 
served spectra of SN 2007if. In the fit, the local continuum 
is estimated separately for each spectrum as a cubic poly- 
nomial, and the absorption is fully modeled using separate 
wavelengths for both Na I D lines using the SNIPS instru- 
mental resolution of 6 A. The resulting probability density 
was then marginalized over the redshift error to obtain the fi- 
nal estimate. Por the Milky Way dust absorption, we derive 
£W(Na I D) = 0.5 l^g j]^ , which under the "shallow" TBC rela- 
tion corresponds to E(B-V) = 0.072, in good agreemen t with 
the reddening from Schlegel. Finkbeiner & Davis! (11998 ). No 
absorption is detected from the host, and we derive an upper 
limit ^^(Na I D) < 0.14 (95% CL), which corresponds to 
E{B-V) < 0.032 even under the "steep" TBC relation. 

A similar discrepancy betw een the Lira re l ation and 
Na I D absorption was noted by lYamanaka et al.\ (120091) for 
SN 2009dc; their work ended up adopting the estimate from 
^^(Na I D). Under the same assumption for SN 2007if, 
host reddening is therefore unlikely to account for a signifi- 
cant share of SN 2007if's deviation from the colors of nor- 



mal SNe la. The color difference is probably real and intrin- 
sic. If the entire color difference were due to reddening with 
Rv ^3.1, this would make SN 2007 if nearly a full magnitude 
brighter at maximum than SN 2003fg. However, given the 
disagreement between the two host reddening estimates, we 
revisit the potential i mpac t of uncontrolled host reddening on 
our results in Section 15.61 The primary analysis to follow in 
our paper assumes zero host galaxy reddening. 

Adopting the magnitudes at maximum light from the 
SALT2 fits to the rest-frame light curve and a CMB-frame 
redshift of zcmb = 0.0731 ± .0013 (including a 300 km s"' 
peculiar velocity) gives a Hubble-flow distance modulus of 
fi = 37.57 ± 0.03 for the co ncordance cos mology {h = 0.71, 
VLm = 0.27,17a = 0.73; see IWrightl [2006h implies absolute 
magnitudes Mb = -20.23, My = -20.39 assuming zero red- 
dening from the host galaxy. This makes SN 2007if the 
brightest super-Chandrasekhar-mass SN la candidate yet dis- 
covered. The SN is at least some 1.3 mag brighter in V 
than a "normal" xi = Q SN la of the same color, and too 
bright for its SALT2 value of xi by about 1.1 m ag (taking 
Mv(xi = 0, c = 0) = -19.07, a = 0. 13, /3 = 1 .77, after lGuv etal\ 
l2007l) . 

The temporal sampling and S/N in ANDICAM J are much 
lower than in the optical bands, due to the limited time avail- 
able in our SMARTS observing program. However, the light 
curve suggests a second maximum starting 15-20 days af- 
ter B-band maximum light, lasting until about 40 days after 
B-band maximum light and then dropping dramatically. The 
height of the first maximum is not clear from the data; the sec- 
ond maximum probably occurs 28-34 days after maximum 
light, at a height of 7 ^ 19. Pigure |6] shows a fit of a fourth- 
order polynomial to the data within the range of phase cover- 
age, showing a possible position for the second maximum. 

The estimate of /-band flux at maximum light is similarly 
problematic. If the 7-band light curve near peak is simi- 
lar to the behavior e.g. of /?-band, extrapolated backwards 
nearly linearly without a pronounced maximum, then maxi- 
mum light might be only as bright as the brightest observed 
point, or /max = 18.4. If the first maximum is more pro- 
nounced, as in the high-^^'Ni light curves of Kasen (2006i), 
then it could be as bright as /max = 17.8. We take these esti- 
mates as representative lower and upper limits. 

Given these estimates of the flux in different bands at B- 
band maximum light, and a spectrum similar to the spectrum 
at +5 days, we estimate a maximum bolometric luminosity of 
(3.22±0.15)x lO'^^^ergs-i. 

3.3. Light curve comparison to SN 2003fg 

Pigure [8] shows a direct comparison between the light 
curves of SN 2007if and SN 2003fg. Since a spectral time se- 
ries which might yield accurate /T-corrections for SN 2003fg 
is not available, we have instead /T-corrected our SN 2007if 
light curve to the observer-frame gri bandpasses of SN 2003fg 
at z = 0.244. 

We can see from this comparison that SN 2007if and 
SN 2003fg have very similar decay behavior SN 2003fg is 
certainly fainter than SN 2007if by about 0.3 mag. The rela- 
tive color and light-curve phase between the two SNe is uncer- 
tain due to the lack of maximum-light coverage in ^-band for 
SN 2003fg. SN 2003fg is probably bluer than SN 2007if, but 
not by more than about 0. 1 mag, or B - V = 0.0 6; this is diffi- 
cult to reconcile with the reported B-V = -0.15 dHowell et al.\ 
i2006; .Hicken et a/... 2007.) . We shall return to the interpreta- 
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Fig. 8. — Direct comparison of SN 2007if multi-band light curves 
with SN 2003fg in the latter's observer frame. Solid circles: 
ANDICAM+QUEST+SNIFS BVR for SN 2007if if-corrected to 
SNLS gri at z = 0.244. Open squares: SNLS gri observer- frame 
(z = 0.244) light curves for SN 2003fg (see Howel l et a/.l20M) . 

tion of possible color differences between the two SNe below; 
it is unclear that the excess should be interpreted in terms of 
dust extinction. 

4. SYNTHESIZED NICKEL MASS IN SN 2()07IF 

The maximum magnitude R= 17.2 observed in the search 
and follow-up images suggest a highly overluminous {My ^ 
-20) explosion for SN 2007if, similar- to SNe 2003fg and 
SN 2009dc. This in turn suggests a very large mass of ^^Ni 
— could it be in excess of Mch? In this section we calculate 
a bolometric light curve of the SN, and from an estimate of 
the bolometric luminosit y at maximum light, derive the ^^Ni 
mass using Arnett's rule (lArnettill982l) . 

4.1. Bolometric light curve 

Bolometric hght curves of SNe la are useful for estimating 
the ^^Ni mass ( lArnettl[T982l) . from the bolometric luminosity 
at peak, and the total ejected mass, from an estimate of the 
energy deposition in the ejecta from decay of radioactive ^^Ni 
and ^''Co dStritzinger et a/.ll2006b . iMaeda & Iwamotol (l2009h 
characterize each SN by the decay time of the bolometric light 
curve. 

To produce a bolometric light curve for a given SN la, 
iStritzinger et al.\ (^06) use a parametrized model to interpo- 
late the multi-band light curve in time, which is then numeri- 
cally integrated in wavelength to find the total flux over bands 
as a function of time. We choose to take a slightly different ap- 
proach using the observed spectra of SN 2007if, similar to that 
used by Howell et al. (2009). For each quasi-simultaneous set 
of ANDICAM BVRI observations, we deredden and dered- 
shift the SNIFS spectrum nearest in time, then multiply it 
by a smooth function (in this case a cubic polynomial) fit- 
ted so that the synthetic photometry from the resulting spec- 
trum matches, in a least-squares sense, the corresponding 



ANDICAM imaging photometry in each band. We then cal- 
culate the bolometric flux as the integral of the SNIFS spec- 
trum over all rest-frame wavelengths from 3100-9000 A. Oc- 
casionally (< 10% of the time) an ANDICAM observation in 
one band is missing due, e.g., to S/N considerations from in- 
termittent clouds or bright moonlight; since the phase cover- 
age is overall still exemplary, such gaps are covered via linear 
interpolation between adjacent light curve points, and a sta- 
tistical error bar assigned to the interpolated value based on 
the S/N of the adjacent points. This never occurs in more than 
one of the bands at a time. As a cross-check on the effects 
of the evolution of spectral features on the bolometric flux 
thus reconstructed, we evaluated each bolometric light curve 
point using the two spectra bracketing it in time, and took the 
difference between the measurements as an estimate of the 
systematic error. This difference was found always to be less, 
and usually much less, than the statistical error on each light 
curve point. 

This handles only the optical section of th e light curve, 
however. Given the discussion in Section \52\ below, we can 
expect much of the flux to be reprocessed into the NIR. To 
account for this, we repeat the above procedure for the NIR, 
by normalizing the synthetic J magnitudes over the spectra of 
SN 1999ee to match the observed /-band photometry. We ex- 
pect systematic errors resulting from this approximation to be 
small, since SN 1999ee is itself an overluminous LVG SN la. 
We then integrate the flux from 9700-24800A. To compen- 
sate for the poor phase coverage and S/N in ANDICAM J, 
we normalize to the fourth-order polynomial fit in Figure |6] 
We estimate, or place limits on, the /-band contribution out- 
side the bounds of phase coverage with a piecewise linear 
extrapolation based on the 7-band light curves from Kasenl 
^006). Most of the high-^''Ni models show a decline of 
0.08 mag day"' after the second peak, which we adopt for 
rest-frame day +55 to day +70. The subsequent behavior is 
uncertain and we allow for the possibility of a decline any- 
where in the range 0.02-0.06 mag day"'. These limits and 
uncertainties are propagated into the bolometric light curve 
error bars. We find that the NIR correction to the bolometric 
luminosity ranges from 30% of the optical luminosity (near 
the second maximum) to 10%-15% after day +60. The NIR 
correction is less than 5% around B-band maximum light, so 
that the expected systematic error on the ^^Ni mass is small. 

The resulting bolometric light curve is shown in Figure |9] 

4.2. Estimate of the ""^Ni mass from Arnett's rule 

To estimate the ""^Ni mass, we calculate the rate of radioac- 
tive energy deposition in the ejecta, 

J^iad = -^56NiA56Ni256NU^ eXp(-A56Njf) 

A56CoA56Ni 



+ A/; 



■ (25<5Co,f+ + 256Co,'y) 



/\56Ni — /\56Co 

X (exp(-A56Cof ) - exp(-A56Nif )) , 



(2) 



where t is the time since explosion, ^6^; =M56is,i/(56 AMU) 
is the number of ^^Ni atoms produced in the explosion, Assn; 
and A56C0 are the decay constants for ^^Ni and ^''Co (lifetimes 



5'>Co,7 



and 



8.8 days and 111.1 days) respectively, and 256^;^, Q 
Q5(,co.e+ are t he energies releas ed in the different stages of the 
decay chain dNadvozhinl 1 9941) . We ado pt an efficiency factor 
a of o rder unity such that Lboi = aLrad dHoflich & Khohklovl 
119961; iJefferv. Branch. & Baron' '2006'). We use a fidu- 
cial value of a = 1.2 (Jefferv. Branch. & BaronI 120061 ; 
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Fig. 9. — Rest-frame bolometric light curve of SN 2007if, synthe- 
sized from SMARTS BVRIJ data, SNIPS spectroscopy, and NIR 
spectraof SN 1999e e taken with ISAAC on the VLT at Cerro Paranal 
dHamuv et a/.l2002h . A fit of the radioactive energy deposition curve 
to the data taken more than 60 days after explosion is shown. 



iHowell et g/.l 120061 \2009h. within the range 0.8 < a < 1.6 
repre sentative of various models dJefferv. Branch. & BaronI 
120061) . Given the quality of our data, the uncertainty on a 
is the limiting systematic error in our analysis. 

Since ROTSE-lllb detected the SN at unfiltered magni- 
tude 19.5 ±0.1 as eai-ly as August 16.3 UT (Akerlof ef aZ] 
120071) . the observer-frame rise time must be at least 20 days. 
Palomar-QUEST detects no object at the SN location on Au- 
gust 9.0, however, placing an upper limit of RG610 > 21.5 
at 95% confidence; thus the observer-fram e rise time cannot 
be any longer than 26 day s. A parabolic (iRiess et g/.l 1 19991 : 
lAldering. Knop. & N ugent 2000) fit to the first four points on 
the ROTSE light curve ( Akerlof ef a/. 2007) suggests an ex- 
plosion date of 2007 August 10.7 (MJD 54322.6); this gives 
a rest-frame rise time of 24.2 ± 0.4 days, which we adopt for 
our analysis. 

Using these numbers we derive a ^^Ni mass of Mssn; = 
( 1 .72 ± 0.09) X (a/ 1 .2) Mq, somewhat larger than SN 2003fg 
and SN 2009dc, and already in excess of the Chandrasekhar 
mass for a = 1.2. This estimate assumes no extinction due 
to dust in the host galaxy. The Chand r asekha r-mass DETl 
model of i Khokhlov. Miiller & HoflicE (12006^. in contrast, 
makes only 0.92 Mq of ^^Ni, approximately the theoretical 
upper limit for Chandrasekhar-mass detonations. We will re- 
turn to the ^''Ni mass estimate below in the context of more 
detailed models of the SN 2007if progenitor. 

5. CONSTRAINTS ON THE SN 2007IF TOTAL MASS AND DENSITY 
STRUCTURE 

There have been several attempts to explain how a 
Chandrasekhar-mass white dwarf in a single-degenerate sce- 
nario could give rise to an unusually luminous explo- 
sion, u^iwny_jel^2IlS_2!L^£E3I?^J£Lfr spherical symme- 
try. iHillebrandt. Sim. & Ropkel (l2007h simulated aspherical 
three-dimensional SN la explosions, which at their brightest 



achieved Mboi = -19.9 at maximum, but at the cost of some 
fine-tuning and requiring an almost complete conversion of 
the star to ^^Ni. They also failed to reprod uce the low ve- 
locities found in SN 2003fg-like SNe. Kas en et aL\ ( 120041) 
considered initially spherical explosions which lost spherical 
symmetry through the interaction of the SN ejecta with a non- 
degenerate companion star, forming a conical hole through 
which deeper, hotter ejecta became visible along certain lines 
of sight; such a model produces a SN with a 1991T-like spec- 
trum, bluer colors and only a modest 0.3 mag) luminosity 
enhancement relative to normal SNe la. Neither of these mod- 
el s are consiste n t with our observations of SN 2007if. 

iTanaka et al.\ (12009 ) present polarimetry observations for 
SN 2009dc which show no significant departures from spher- 
ical symmetry. This would seem to imply that, whether or not 
the progenitor star is aspherical, a highly aspherical explo- 
sion is not necessary to attain a very high luminosity. How- 
ever, it is reasonable to ask whether observations can pro- 
vide further constraints on the mass of the progenitors in such 
systems apart from measuring the SN's luminosity and ve- 
locity near peak, within the context of an explosion with no 
large-angular-scale anisotropics. One way of approaching this 
question is via the late-time bolometric light curve (Jeffe^ 
1999; Stritzinger et al. 2006), which requires an accurate es- 
timate of the kinetic energy of the explosion and the corre- 
sponding velocity scale for the ejecta density profile. 

Our measurements (see Figure |4| show an unusually low 
Si II A6355 velocity in the ejecta, as has been seen in other 
super-Chandrasekhar SN la candidates. We also see a fairly 
flat velocity evolution, consistent with a plateau in the Si II ve- 
locity (Sv < 500 km s"'; lOuimbv. Hoflich & Wheeleiil2007h 
starting as early as day -9 and extending at least as late as day 
+7. These low velocities have previous ly been attributed t o 
large binding energies for these events (Ho well et al.\ \2006). 
In the following section we explore an alternative — that the 
low velocity might result from deceleration of the outer layers 
of ejecta by a massive envelope surrounding the SN progen- 
itor, producing an overdense shell in the ejecta. We present 
a simple parametrized treatment of the shell structure which 
allows us to constrain the true underlying kinetic energy of 
the explosion, the mass of the shell and the envelope, and the 
total ejected mass. 

5.1. Shell structure in SN 2007if 

lOuimbv. Hoflich & Wheeleii (I2007h studied the influence 
of shell structure in SN la ejecta for a normal SN la event, 
SN 2005hj. They pointed out that an overdense spherical shell 
in the ejecta expanding at a given velocity slows the reces- 
sion of the photosphere, producing a plateau at that veloc- 
ity in the time evolution of the widths of absorption features, 
e.g.. Si II. The shell would also be characterized by a duration 
and a velocity-space width from the spectra. Photometrically, 
SNe la with an overdense shell in the outer layers should ap- 
pear redder and brighter than SNe la with no shell. Existing 
models of shell-structure SNe la have been calculated only 
for Chandrasekhar-mass explosions, but one should expect 
analogous physical effects in super-Chandrasekhar-mass ex- 
plosions. Such models provide an alternate interpretation for 
the low Si II velocities seen in SN 2003fg and SN 2007if, al- 
though a large progenitor binding energy may also contribute. 

Shell structure occurs naturally in DD merger scenarios, in 
which the tidal disruption of one or both white dwarfs creates 
a common envelope of carbon and oxygen around the cen- 
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tral merger product. Collision of the ejecta with the envelope 
then accelerates the envelope to higher velocities and creates 
a strong reverse shock which creates a shell in the outermost 
SN ejecta. Hoflich & Khohklov (1996) simulate several in- 
stances of common envelopes in their DET2ENVN "tamped 
detonation" (TD) models, which place two merging white 
dwarfs of total mass 1 .2 solar masses in an envelope of mass 
OAxNMq. DET2EN V4 and DET2EN V6 produce B - V - 
0.2 near maximum light, significantly redder than SN 2006gz 
but consistent with 2007if. The DET2ENV6 model, with a to- 
tal system mass of 1.8 Mq, also has AniisiV) = 0.63^qq9, the 
slowest of any model in that work, compared with Amig(y) = 
0.50±0.07 for SN 2007if. (The numberfor SN2007if is inter- 
polated from the observations using the SALT2 fit; the domi- 
nant source of error for both the models and the observations 
comes from the uncertain phase of maximum light). 

The single-degenerate pulsating delayed detonation (PDD) 
models of Hoflich & Khohklov (1996) also have shell struc- 
ture in the ejecta. Here an initial deflagration phase pre- 
expands the white dwarf; the outermost layers begin to rec- 
ollapse and are then entrained into a shell by the underly- 
ing ejecta from a subsequent detonation. These models do 
exhibit an intrinsically red B-V color at maximum, though 
with a wide spread (0.05-0.60). However, they have unusu- 
ally short rise times (13-15 days) and slower dechnes than 
observ ed for SN 2007if. Most of the iHoflich & Khohklovl 
( 119961) models have rise times which are too short — around 
15 days in V, compared to the observed typical value of 
19 days ("Conley et flZ."2006). This aspect of the models could 
no doubt be improved for a more accurate match to obser- 
vations. Nevertheless, it is striking that the shell-structure 
explosion models reach the extremes of the rise time distri- 
bution with fairly small dispersion, with most PDD models at 
13 days and the TD models as large as 22 days. The one PDD 
model that achieves a sufficiently long rise time and slow 
decline, PDD535, is much fainter even than normal SNe la 
(My = -17.77) and has an extremely red color of B-V = 0.60. 
This model produces about 20% of the kinetic energy of the 
DET2ENVN models. 

A separate question arises regarding the widespread use 
of the differentially rotating models of lYoon & Langed 
(2005) in treatments of super-Chandrasekhar-mass S Ne la 
jHowefl et all 120061; iJefferv. Branch. & BaronI l2006t 
iMaeda & Iwamoto' ""20091), including this paper (see be- 
low). Piro (2008) found that inclusion of the baroclinic 
instability and of magnetohydrodynamic processes in white 
dwarf interiors inhibits the Kelvin-Helmholtz instabilities 
expected in differentially rotating white dwarfs, limiting 
such white dwarfs to rigid rotation and apparently making 
single-degenerate super-Chandrasekhar-mass white dwarfs 
untenable. Even discounting such effects, however, su- 
permassive white dwarfs accreting from a non-degenerate 
compa nion are not exp ected to evolve to masses above 1.7 
Mq ( Chen & Lil 120091) ; the progenitor of SN 2007if was 
almost certainly more massive than this. 

The interpretation of a velocity plateau as the sign of an in- 
teraction with an envelope containing a few tenths of a solar 
mass of material also poses problems for models involving the 
explosion of single-degenerate supermassive white dwarfs. 
To the extent that such an envelope is necessary to explain the 
kinematics of 2003fg-like SNe la, it must be composed of car- 
bon and oxygen; a hydrogen/helium envelope of comparable 
mass would most likely produce strong emission lines such 
as those seen in the "SNe 11a" 2002ic (.Hamuv et ai. .2003; 



|Wood- Vas ey TFoZl l2004h and 2005gj dAldering et al.\ 120061; 
Prieto et a/. l 2007h . 

The above considerations disfavor existing single- 
degenerate explosion models, although it is conceivable that 
some PDD model could be tuned to fit our observations. The 
TD models provide a better fit to the observations. For further 
analysis below, we therefore discount models of pulsating 
delayed detonations for the time being, and concentrate in- 
stead on tamped detonations representing double-degenerate 
mergers, using DET2ENVN as our main point of reference. 

5.2. Near-IR light curves and mixing of the ejecta 

Our estimate of the tota l ejected mass will follow ' Jeffervl 
(1999) and Stritz inger et a l. (2006), relying on the estimated 
transparency of the ejecta to gamma rays from ^''Co decay. 
The model requires an assumption about the distribution of 
^^Ni in the ejecta. Although this incurs some uncertainty, we 
argue that the ejecta of SN 2007if are fairly well mixed based 
on the near-infr ared (NIR) light curves. 

Kasen (2006) presents radiative transfer Monte Carlo cal- 
culations which illustrate the relevance of mixing in SN la 
ejecta to the appearance of an inflection point or second max- 
imum in their NIR light curves. The recombination of iron 
provides an efficient means of redistributing light from UV 
and blue bands to the NIR, and recombination happens at a 
more or less sharp front in SN la ejecta as they expand and 
cool. Thus in SNe with compact cores of iron-peak elements, 
a wave of recombination of Fe III to Fe II powers the first NIR 
light curve peak, while the recombination of Fe II to Fe I pow- 
ers the well-known second peak. Supernovae in which the 
ejecta are well-mixed have less pronounced second maxima 
in the near infrared than those in which the ejecta are strati- 
fied. Supernovae which produce less ^^Ni also show less pro- 
nounced second maxima, and these maxima occur earlier in 
phase and are more closely spaced with respect to the first 
maximum, since the temperature of the ejecta is lower and 
recombination of iron-peak elements occurs earlier. 

The /-band light curve of 2007if shows only a very slight in- 
flection which is difficult to measure. The 7-band light curve 
shows stronger, but still quantitatively weak, evidence for a 
second maximum. In the context of this model, the observed 
/- and 7-band light curves of SN 2007if suggest either a low 
^^Ni mass, ejecta which are highly mixed, or both. Assuming 
stratified ejecta, the shape, timing and duration of the /-band 
second peak is broadly consistent with a large amount of ''^Ni 
(0.6 Mq or more). Given the luminosity of the event, a plau- 
sible interpretation is that the ejecta contain a large amount of 
^^Ni, much of which lies near the surface — either because 
it has been mixed to higher velocities in the explosion, or be- 
cause the ejecta are composed mostly of ""^Ni and little else. 
This interpretation is also supported by the faster decline in 
Z?-band relative to SN 2003fg, suggesting that the Z?-band flux 
is indeed being efficiently reprocessed to NIR wavelengths 
dKasen & WooslevI l2007h ; a smaller amount of ^^Ni might 
manifest as a sharper decline in all bands (i.e. the bolometric 
light curve) rather than B-band only. It is also consistent with 
the appearance of Fe III lines near maximum, as mentioned in 
SectionO 

Khokhlov, Muller & Hoflich (120061) point out that in SNe 
with shell structure, the shock interface between the shell 
and the high-velocity ejecta is Rayleigh-Taylor unstable. This 
could provide a mechanism for mixing of iron-peak elements 
into the shell in SN 2007if. If so, we might expect the asso- 
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ciated features — a weak NIR second maximum and a rel- 
atively sharp B-band decline rate — to be a generic feature 
even of shell-structure SNe with less massive progenitors. We 
might therefore also expect the usual width-luminosity rela- 
tionship to break down for shell-structure SNe la, produc- 
ing SNe which are systematically brighter for their (B-band) 
declin e rate than SNe la with near-exponential density struc- 
tures ( Kasen & Wooslevll2007l) . Whether the bolometric light 
curves still followed a width-luminosity relationship would 
depend in part on the efficiency of the optical-NlR reprocess- 
ing. 

5.3. Mass estimates and photospheric evolution for SN la 
models with exponential density profiles 

In the optically thin limit, and when the gamma rays come 
from a single radioactive species. Equation ^ is modified by 
multiplying the Qco.^ term by a factor 1 - exp(-T^), where 
is the optical depth to scattering of gamma rays from the decay 
of ^*Co. Such approximations are useful in des cribing the 
state of the eject a by day +60 after explosion (Jeffervl l 19991: 
IStritzinger et a/.ii200 6). by which time less than 0.5% of the 
initial ^^Ni should remain. 

The optical depth should scale as f"^ for SN ejecta in 
homologous expansion. The evolution of the ejecta density in 
a SN la can be well described by 

P(v,f) = Pc,()(fo/f)^exp(-y/y,), (3) 

where pc,o is the central density at fiducial time fo and is the 
scaling velocity. By taking fo to be the time at which = 1, 
i.e. the transition to the optically thin regime (Jefferv 1999) 
and integrating outwards from the center of the SN, we ob- 
tain pc,o = (K7^Vcfo)~\ where k,^ is the opacity for gamma-ray 
absorption, and <7 is a form factor describing how the ^^Ni is 
distributed throughout the ejecta. The total mass is 

MwD = 87^pc.o(v^;^o)^ 
and the kinetic energy is 

EK = 6MyNDv]. 

We can then use these relations to estimate the total mass as 

MwD=— (v,fo). (4) 

In their analysis of a sample of normal SN la bolometric light 
curves, Stritzing er et a l. (2006) adopt q = 0.33 ±0.09 (appro- 
priate for ejecta in which ^^Ni is mixed evenly throughout), 
K,^ = 0.025 ±0.03 g cm"^, and = 3000 km s~'. They empha- 
size that the procedure tends to produce low ejected masses 
(< 1 .0 Mq) for many SNe la, and comment on the accuracy of 
the input parameters. The values of and q are not expected 
to vary much, and the largest single uncertainty is expected to 
arise from the value of v^. used, which may have uncertainties 
on the order of 20% and is squared in the above expression. 

For ejecta following an exponential density distribution, the 
photospheric velocity Vphot can also be calculated given v,. by 
solving Topt(vphot) = 2/3 (~ 50% escape probability). This 
leads to 

/3KoptMwD\ 

where Kopt is now the opacity for optical photons. For 
a near-exponential density profile, we would expect Si II 



A6355 to track the photosphere (and hence v^) for observa- 
tions near maximu m light; this is the assumption made by 
iMaeda & Iwamotol ( l2009i) in their parameter study of super- 
Chandrasekhar-mass SNe la. 

5.4. Incorporating the effects of a shell 

The ejected mass estimate in Section |53] hinges on the as- 
sumption that the ejecta follow an exponential density dis- 
tribution with velocity. This is no longer true for SNe with 
shell structure. However, the DET2ENVN models have an es- 
sentially undisturbed exponential density structure in the cen- 
trally concentrated ejecta underneath the shell. Although the 
presence of a shell modifies the velocity structure and colors 
of the SN considerably, it should have negligible effect on 
the bolometric light curve at sufficiently late times, since the 
opacity is an order of magnitude smaller than ^opt in most 
realistic scenarios. Due to geometrical dilution, the shell be- 
comes transparent to gamma rays by maximum light, after 
which the exponential models should apply. The shell mass 
should consist simply of the ejected mass above some velocity 
v'sh, which looks like 

MwD V J 

where Q{a,x) = 7(a,x)/r(a) is the incomplete gamma func- 
tion. Similar gamma-function integrals can be calculated to 
give the fractions of the optical depth, momentum, and ki- 
netic energy, for which the values of a are 1, 4, and 5, respec- 
tively. For a typical model with Ve = 2750 km s"', the material 
above 9000 km s"' carries about 35% of the ejected mass and 
75% of the kinetic energy, yet contributes less than 5% to the 
gamma-ray optical depth as measured from the center of the 
explosion to infinity. 

We therefore expect that the bolometric light curve will 
yield a reasonable mass estimate for SN 2007if, provided 
we choose Vg appropriately. However, in the presence of a 
shell, the Si II velocity is not representative of the kinetic en- 
ergy velocity Vg during the velocity plateau, but instead rep- 
resents the shell velocity. Our observations show no signifi- 
cant change in the Si II velocity over a 16-day period (-9 days 
to +7 days), inconsistent with the photospheric evolution of 
a near-exponential density profile even for SNe la with very 
low kinetic energy. Although Si II A6355 is sometimes used 
at late times (+40 days) to estimate v^, the line is too heav- 
ily blended with other features in the SN 2007if post-plateau 
spectra to measure its velocity. We therefore cannot observe 
v'e directly in our work here, and instead marginalize over it 
as a nuisance parameter. 

For double-degenerate progenitor models with envelopes, 
we can derive the shell and envelope masses given via con- 
servation of momentum and optical transparency arguments. 
Shortly after explosion, all SN ejecta above v = v'sh are coast- 
ing freely and the envelope with mass Mgnv is stationary. At 
late times, the envelope has been accelerated to higher ve- 
locities, and the outermost SN ejecta are moving together in 
a shell of negligible thickness with velocity Vsh (which will 
match our measured Si II velocities). In the DET2ENVN 
models, the specific momentum in the C/O envelope is some 
50% higher than in the shell. Conservation of momentum then 
gives 

Pinitial = 3MwdVc2 ( 4, — ) = Pfinal 
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MwdG 3, 



+ 1.5Me, 



Vsh- 



Some algebra then gives 



Me, 



Vsh \ Ve J \ Ve 



WD- 



(7) 



(8) 



The sum Mtot = Mavi+My^D should then equal the mass of the 
two white dwarfs in the original close binary system. The 
combined envelope-plus-ejecta shell should remain opaque 
for a time 

u. = ^^^^. (9) 



Vsh 



We can calculate fsh, assuming fractions /ipe and /ime of iron- 
peak (Ni, Co, Fe) and intermediate-mass elements (including 
carbon and oxygen) respe ctively, by using the approximate 
line-opacity prescription of lMazzali et al\ (1200 lb . 



/vopt = 0.5 [0.25/ipe + 0.025/ime] cm^ g"' (10) 

(see also ' Maeda & Iwamotol l2009). Our use of this relation 
assumes the material in the shell has been mixed. The con- 
served momentum and the observed transparency of the shell 
provide indirect constraints on the underlying, unobserved Vg. 

5.5. Confidence regions on the progenitor mass 
In our analysis, following iJeffervl (Il999h . we adopt = 
0.029 cm^ g~' and q = 0.33 (since an even distribution of 
^^Ni is supported by our observed NIR light curves). The 
appropriate value for will depend on the kinetic energy, 
which in turn will depend on the mass and on the composi- 
tion of the ejecta. To get a self-consistent mass estimate, we 
generate a suite of semi-analytic progenitor models similar 
to those presented in Jefferv, Branch. & BaronI (12006), but in- 
cluding an envelope with characteristics similar to those found 
in DET2ENVN. We parametrize our models by the white 
dwarf mass Mwd and the fractions f, of ^''Ni, other iron-peak 
elements ("Fe"), intermediate-mass elements ("Si") and un- 
burned carbon and oxygen ("C/O"), the radioactive energy 
deposition parameter a, and the ratio /e„v of the C/O enve- 
lope mass Menv to the ejected mass Mwd which contributes 
to the original binding energy. The binding energy was cal- 
culated from the differ entially rotating white dwarf models of 
lYoon & Langeii (12005 ). and the kinetic energy was calculated 
as the difference between the binding energy and the nuclear 
energy released. 

After calculating for each model, we could then make 
predictions for the observed bolometric light curve to be com- 
pared directly with observations. In matching the bolometric 
luminosity at peak (Arnett's rule) jointly with the other ob- 
servations, we sample values of a between 0.8 and 1.6 to 
cover a reasonable range of potential variation in different 
explosion scenarios. However, for our primary analysis we 
adopt a fiducial value of a = 1.3, characteristic of the tamped- 
detonation models DET2ENVN w hich lie within ±0.05 of 
this value ( Hoflich & Khohklov '1996). We calculate Vsh given 
Vf and Menv by solving Equation ([8j numerically, and we cal- 
culate fsh from equations |9] and [TO] (using /ipe = /fb + /ni and 
/ime = /si + /co)- In constraining the models we require that 
Vsh match the observations (9000 ± 500 km s"' ) and that fsh be 
longer than the observed lifetime of the plateau (+7 days af- 
ter B-band maximum, or +32 days after explosion). Although 



the spectral characteristics begin to change at day +10 after B- 
band maximum, including a weakening of Si II A6355 which 
signals the end of the plateau phase, we do not enforce a hard 
upper limit on the length of the plateau. 

The value of Mtot was sampled uniformly between 1 .4 and 
2.8 Mq, with the highest allowed mass corresponding to 
the merger of two Chandrasekhar-mass white dwarfs. We 
chose 2 Mch as the most conservative upper bound on the 
system mass, since observations, or the theory of common- 
envelope evolution thought to produce double-degenerate sys- 
tems, do not yet provide_secure constraints (see for 
ple iDobbie et a/.l 120061 [200^ . The values /pe, /m, /si, and 
/co (representing the composition of the bound progenitor 
mass) were sampled uniformly between and 1 such that 
their sum was unity. The value of /env = Mgnv/MwD was 
varied independently between and 0.5. The central den- 
sity of the merged progenitor was fixed at 3 x 10'' g cm"-*, 
higher than the DET2ENVN models but typi cal of central 
densi t ies of white dwarfs assum ed elsewhere dHowell et al.\ 
l2006t iMaeda & Iwa"rnotoll200Si) . For models with /env > 0, 
i.e., models with shells, we assume that all of the IME are 
mixed into the shell, to agree with our observations of very 
weak or absent Si II in post-plateau spectra. 

To generate confidence regions, we evaluated a chi-square 
(X^) for each progenitor model when compared to all avail- 
able observations. Each observation — the maximum bolo- 
metric luminosity, bolometric luminosity measurements more 
than 60 days after maximum, and velocity measurements dur- 
ing the plateau phase — contributed equally to according 
to its error. Each model was then weighted by the probability 
of observing the calculated value of x^, and binned in a his- 
togram. The confidence regions drawn contain the stated frac- 
tions of the total probability of all models calculated. There 
may be additional constraints on the actual distribution of 
these parameters realized in nature, which we do not include 
in the sampling or weighting of different possible combina- 
tions of parameters. For example, SN la explosion physics 
should constrain the relative abundances of various elements 
in the ejecta; our constraints rely only on the kinematics of 
the explosion given the nuclear energy released. Due to this 
and other potential limitations, we therefore do not claim to 
pick out a single set of values which is most likely, nor to re- 
produce the detailed shapes of the distributions of underlying 
quantities. However, this approximate method should allow 
us to rule out untenable models. The size of the confidence 
regions, interpreted loosely, should reflect the uncertainty on 
the measurements themselves. 

Figure [To] shows the first-pass results of this procedure for 
the case a= 1.3. The constraints are surprisingly tight: the 
data are consistent with a model containing 2. 18 ±0. 15 Mq of 
central ejected material, with 1 .55 ± 0.09 Mq of the mass be- 
ing ^''Ni, wrapped in a C/O envelope of mass at least 0.3 Mq. 
The shell of reverse-shocked SN ejecta contains about a third 
of the central ejected mass, or some 0.7-0.9 Mq of material. 
The minimum chi-square value is xt =25.1/21 = 0.93, indi- 
cating that the model is able to satisfy all available constraints. 
Chandrasekhar-mass progenitors, and progenitors with no en- 
velopes, are ruled out at high significance. We place a lower 
limit of 2.05 Mq (99% CL) on the total mass of the system, 
and 0. 1 Mq at similar confidence on the mass of the envelope. 

The upper-right-hand plot in Figure [10] illustrate the trade- 
off between unburnt, initially bound carbon/oxygen in the 
SN ejecta and more loosely bound carbon/oxygen in the sur- 
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TABLE 3 

Variation in constraints with model-dependent parameters 



a 




Rv,host Pc (lO' gem 3) 


M,o, (Mq) 


MwD (Mq) 


Mmi (Mq) 


/env 


/sh 


Prob." 


0.80 


0.00 


— 3.00 












< 10-5 


0.90 


0.00 


— 3.00 


2.26 ±0.06 


2.00 ±0.06 


1.84 ±0.06 


0.14±0.01 


0.39 ±0.02 


0.004 


1.00 


0.00 


— 3.00 


2.28 ±0.09 


2.03 ±0.08 


1.78 ±0.07 


0.14±0.02 


0.38 ±0.02 


0.078 


1.10 


0.00 


— 3.00 


2.32±0.12 


2.08 ±0.11 


1.71 ±0.09 


0.15±0.03 


0.37 ±0.02 


0.378 


1.20 


0.00 


— 3.00 


2.36±0.15 


2.13±0.13 


1.63 ±0.09 


0.15±0.03 


0.36 ±0.03 


0.798 


1.30 


0.00 


— 3.00 


2.41 ±0.17 


2.18±0.15 


1.55 ±0.09 


0.16±0.04 


0.34 ±0.03 


1.000 


1.40 


0.00 


— 3.00 


2.45 ±0.18 


2.23 ±0.15 


1.48 ±0.08 


0.16±0.04 


0.32 ±0.03 


0.836 


1.50 


0.00 


— 3.00 


2.49 ±0.1 8 


2.28±0.16 


1.41 ±0.08 


0.17±0.05 


0.31 ±0.03 


0.488 


1.60 


0.00 


— 3.00 


2.53 ±0.17 


2.33±0.15 


1.34±0.07 


0.18±0.05 


0.30 ±0.03 


0.208 



1.30 


0.09 


2.1 


3.00 


2.42 ±0.11 


1.30 


0.09 


3.1 


3.00 


2.45 ±0.09 


1.30 


0.18 


2.1 


3.00 


2.50 ±0.06 


1.30 


0.18 


3.1 


3.00 





2. 17 ±0.09 1.87 ±0.09 0.15 ±0.02 0.36 ±0.03 0.546 

2.20±0.08 1.98±0.10 0.15±0.02 0.36±0.03 0.231 

2.25 ±0.06 2.03 ±0.10 0.14 ±0.02 0.36 ±0.03 0.019 

— — — < 10-5 



1.30 0.00 — 0.03 2.72 ±0.07 2.40 ±0.06 1.61 ±0.08 0.39 ±0.03 0.40 ±0.03 0.213 

1.30 0.00 — 0.10 2.66±0.11 2.36±0.09 1.61 ±0.09 0.31 ±0.04 0.39±0.03 0.631 

1.30 0.00 — 0.30 2.62±0.14 2.33±0.11 1.60±0.09 0.25±0.04 0.38±0.03 1.044 

1.30 0.00 — 1.00 2.55±0.16 2.29±0.13 1.58±0.09 0.20±0.04 0.37±0.03 1.332 

Note. — Constraints on properties of models as a function of the (model-dependent) pai'ameters a =Lbol/^'rad: £(B~^)host> Rvhost, 
and pc- Columns, left to right: a; assumed host extinction E(B-V) and Ry; assumed central density pc of the white dwarf before 
explosion; mass of the white dwarf merger product Mwd; '^Ni mass synthesized in the explosion; ratio of envelope mass to white 
dwarf mass; ratio of shell mass to white dwarf mass; total probability, marginalized over all free parameters. 

" Normalized to a = 1.3, no host extinction, pc = 3 X lO' g cm"' (fiducial analysis). 



rounding envelope. An envelope turns out to be much more 
effective than bound unburnt material in producing the low 
velocities observed in the SN. The slope of the contour shows 
that about 4 times as much bound material is required to re- 
main unburnt (decreasing the released nuclear energy for a 
given binding energy) as is needed to produce the same effect 
in an envelope (decelerating initially fast-moving ejecta in a 
collision). The inferred shell/envelope mass is therefore rel- 
atively insensitive to the amount of material left unburned in 
the supernova explosion. 

Given that C II A6580 is the only visible carbon line in 
the velocity-plateau spectra around Z?-band maximum, and 
is detected only weakly, it seems unlikely that a great deal 
of unburned carbon/oxygen is mixed into the shell. Our ob- 
servations can probably be adequately reproduced by setting 
fco = 0, so that all the unburned material is required to lie 
in the envelope (which must then have a mass of around 0.3 
Mq). This is more in line with the hydrodynamic behavior of 
the DET2ENVN models, and with our understanding of the 
physics of normal SNe la, in which almost all of the material 
in the original white dwarf is burned in the explosion. Fig- 
ure [TT] shows how the mass constraints change when we set 
fco = 0; the low-^^Ni models composed mainly of nickel and 
carbon/oxygen are excluded here. When the whole of Mwd 
must be burnt, much more additional mass must be included 
to increase the binding energy sufficiently to reproduce the 
appropriate value of Vg. Under these assumptions, the 99% 
CL lower limit on the system mass shifts upward to 2.20 Mq. 

The results of varying a are shown in Table [3] The allowed 
region of parameter space shows a tight anticorrelation be- 
tween the central white dwarf mass and the ^''Ni mass. This 
intriguing behavior appears to result from the shape of the 
bolometric light curve: models with lower ^^Ni mass favor 
higher values of to, i.e., the gamma-ray escape fraction at a 
given time is lower in these models. It is interesting that the 
data seem to prefer a = 1.3 even with no outside input, pro- 
viding further confirmation that this is a good choice for in- 



terpretation of the data. 

5.6. The impact of uncertain reddening 

The analysis in the previous section was carried out assum- 
ing E(B - y )host = 0. If the SN is reddened by dust, we might 
expect the luminosity and the ^^Ni mass to be higher, and to 
produce a higher total mass. The other observations limit the 
extent of this effect, since an increased mass also leads to 
changes in the kinetic energy and density structure. Given 
the hard upper mass limit of 2 Mch in our modeling, apply- 
ing an uncertain reddening correction has little impact on our 
final mass constraints, but if certain reddening scenarios are 
shown to be incompatible with our observations, we can rule 
them out, even without recourse to the upper limit on Na I D 
absorption. 

To constrain possible reddening scenarios, we repeated the 
analysis of Section |53] for the case a= 1.3, generating new 
bolometric light curves which had been reddened according 
to various prescriptions. The results are shown in Figure [T2I 
These contours represent, in effect, different slices through 
the probability density in parameter space. We normalize the 
probability density displayed in each contour plot to the total 
probability in the zero-reddening case. This allows us to see 
how the relative total number of grid points which satisfy the 
observations changes as the reddening is varied. For example, 
the 68% contour in each of the plots is set at a level which 
encloses 68% of the probability in the zero-reddening case. 

We find that increasing degrees of reddening shift the al- 
lowed contours to higher ^^Ni mass, as expected, although 
the total mass is not strongly affected. Our observations indi- 
cate Av < 0.38 (99% CL), and are therefore incompatible with 
cases in which the Lira relation holds with E(B-V)ho>it = 0.18 
and either Ry = 2.1 01 Ry = 3.1. High extinction tends to re- 
strict solutions to a very small region of parameter space, in 
which the ejecta are composed almost entirely of ^^Ni. In- 
termediate cases with E(B-V) = 0.09 are consistent with the 
data and possibly also marginally consistent with the Lira re- 
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Fig. 10. — Constraints on progenitor models for SN 2007if, assuming a = 1.3 as for DET2ENVN and allowing for a floating fraction of 
unbumed carbon/oxygen mixed into the shell of reverse-shock material. Colored regions are 68% (red), 90%, 95%, 99%, and 99.7% CL (blue). 
Top left: white dwarf mass Mwd vs. nickel mass Mse^i in solar masses. Top right: envelope mass fraction /e„v vs. carbon-oxygen fraction /co- 
Bottom left: probability distribution for the total mass Mmt = Mwd +Mcnv. Bottom right: probability distribution for the fraction of ejected mass 
/sh = Msh/MwD of material in the dense shell of reverse-shocked ejecta. 



lation at the extremes of systematic deviation within the set 
for which it was validated (Folatelli et al. 2009). 

We conclude that, while some of the red color of SN 2007if 
may be put down to dust, it is likely that the Lira relation over- 
estimates the extinction, although the slope of the color evo- 
lution of SN 2007if is similar to normal SNe la in the range of 
light curve phases where the Lira relation holds. If true, this 
has impact for the other observed super-Chandrasehkar-mass 
SN la candidates in which the Lira relation is used to correct 
for reddening, producing large uncertainties in the mass esti- 
mate from the SN luminosity alone for these SNe. 

5.7. The impact of uncertain central density 

One can also vary the central density pc of the SN progen- 
itor at ignition, another input which we held fixed in our cal- 
culations above. The central density for carbon ignition in a 
Chandrasekhar-mass white dwarf is expected to be around 3 x 
10*^ g cm-3 dJefferv. Branch. & Baronll2006h. and is our fidu- 
cial value in the above analvsis. iMaeda & Iwamotol (120091) 
examine central densities as high as lO'" g cm"^, above 
which electron capture is expected t o make collapse to a neu- 
tron star mo re likely than a SN la dNomoto & Kondolll99l !: 
lYoon & Lan ger 2005). The DET2ENVN models, which form 
the primary point of reference for our analysis in the context 
of a tamped detonation, actually have a much smaller central 



density, around 4 x 10^ g cm"^. 

There is therefore a dynamic range of at least two and a 
half orders of magnitude within which pc could be expected 
to vary. As with reddening, decreasing the central density 
can only increase the total reconstructed mass (by decreasing 
the specific binding energy), and so we expect it to have little 
impact on our final mass constraints. However, by seeing how 
our observations constrain allowed masses in different central 
density scenarios, we may be able to place a lower limit on 
the allowed density. 

To this end, we repeat our analysis for£'(B-y) = 0, a = 1.3, 
sampling a range of central densities from 3 x 10^ g cm"-^ 
(close to that assumed for DET2ENVN) to our fiducial den- 
sity. We normalize the probability densities to our fiducial 
model as in Section 15761 above. The results are shown in Fig- 
ure [13] We find that although very low values of the central 
density typical of DET2ENVN are disfavored, a wide range 
of possible values are compatible with the data. While the 
""^Ni mass remains unaffected, the total mass increases, and 
the envelope mass fraction needed to decelerate the explosion 
to the observed velocities increases dramatically. At central 
densities typical of DET2ENVN, an envelope of nearly 1.0 
Mq is needed to match the observations. 

It is possible that other observables associated with our 
spectral time series (for example, the detailed chemical com- 
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Fig. 1 1 . — Constraints on progenitor models for SN 2007if, assuming q = 1.3 as for DET2ENVN and assuming tiiat all carbon/oxygen must 
reside in the loosely bound envelope. Colored regions are 68% (red), 90%, 95%, 99%, and 99.7% CL (blue). Top left: white dwarf mass 
MwD vs. nickel mass MssNj in solar masses. Top right: probability distribution for the envelope mass fraction /env Bottom left: probability 
distribution for the total mass Mmt = Mwd +Me„v. Bottom right: probability distribution for the fraction of ejected mass /si, = Msh/Mwo of 
material in the dense shell of reverse-shocked ejecta. 



position of the ejecta) could provide additional constraints on 
the central density. Such detailed modeling is beyond the 
scope of this paper, and awaits a refined theoretical descrip- 
tion of super-Chandrasekhar-mass SN la progenitors, includ- 
ing hydrodynamic and radiative transfer simulations of the re- 
sulting explosions. 

6. ARE SHELLS COMMON IN SUPER-CHANDRASEKHAR-MASS 
EXPLOSIONS? 

Our detailed observations of SN 2007if have produced a 
relatively specific physical picture of the explosion, providing 
a framework within which to interpret SN 2003fg-like SNe as 
a subclass. A summary of the observational properties of the 
various super-Chandrasekhar-mass SN la candidates found to 
date is given in TableHl we discuss what is known about these 
SNe below, in the context of our findings. The data have been 
put on a uniform distance scale, and measurements and uncer- 
tainties for some quantities have been improved or corrected 
relative to the original sources employing the techniques we 
have applied to SN 2007if. 

Maeda & Iwamoto (20091) argue that an aspherical explo- 
sion is needed to explain the high luminosity of SN 2003fg, 
based on the relatively fast decay time of its bolometric 
light curve (f+i/2 =13 days) and low peak velocity, and raise 
doubts about the high mass of the progenitor as interpreted in 
the binding-energy framework of .Howell et ai. (2006.) . Con- 



versely, they argue that the blue observed color, high veloci- 
ties, and long bolometric decay time of SN 2006gz are consis- 
tent with a roughly spherical super-Chandrasekhar-mass pro- 
genitor. 

Our above a nalysis of SN 2007if all ows us to address 
the criticism of ' Maeda & Iwamotol (12009 ^ regarding the fail- 
ures of a spherical geometry to describe the properties of 
SN 2003fg-like SNe. Their work assumed an essentially 
undisturbed ejecta structure close to an exponential in veloc- 
ity, in which the peak Si II velocity was a suitable proxy for 
the kinetic energy of the explosion. We have seen that, for 
density structures with shells, this assumption is no longer 
true; powerful explosions with high kinetic energy can have 
low-velocity, slowly-evolving photospheres. 

A similar scenario may hold for SN 2003fg. Although it 
is difficult to test because of the limitations of the data, the 
existing observations of SN 2003fg are broadly consistent 
with a SN 2007if-like scenario. The single published spec- 
trum shows a low velocity which may or may not have been 
part of a plateau. SALT2 fits to various near-maximum sub- 
ranges of the SNLS photometry result in values of B-V at 
Z?-band maximum ranging between 0.03 and 0.14, at odds 
with the published estimate B-V = -0.\5 (iHowell etal\200(k 
iHicken et aZ.ll2007h . How ever, the publi s hed (s tatistical only) 
error bars on B-V from iHowell et al.\ (120061) probably un- 
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Fig. 12. — Constraints on progenitor model masses for SN 2007if assuming a = 1.3 (see e.g. upper left panel of Figure [TO]l for different 
possible reddening scenarios. Probability density is normalized to the zero-reddening case. Colored regions are 68% (red), 90%, 95%, 99%, 
and 99.7% CL (blue). Top left: E(B-V) = 0.09, Rv = 2.1. Top right: E(B-V) = 0.09, Rv = 3.1. Bottom left: E(B-V) = 0.1S, Rv = 2.1. Bottom 
right: £(S-V) = 0.18, Rv = 3A. 



derestimate the true error, given that there was no rest-frame 
B-band data at maximum light for SN 2003fg (Howell 2009, 
private communication). A SALT2 fit to all available SNLS 
data gives c = 0.22 ± 0.03, or B - V = +0. 14 at B-band maxi- 
mum after correction for Galactic extinction; this color is con- 
sistent with the color we obtain from synthetic photometry of 
the pu blished SN 2007if specti'um from SNLS (" Howell et all 
120061) . While the evidence is far from conclusive, the spec- 
tral match to SN 20 07if near maximum light, and the argu- 
ments presented in Maeda & Iwamotg (120091 ). suggest that 
SN 2003fg may be best interpreted as a super-Chandrasekhar- 
mass SN la with dense shell structure, like SN 2007if. 

Likewise, SN 2009 dc was also a very re d SN with slow ve- 
locity evolution. Yamanaka e t al.\ (120091) note that the very 
large value of the reddening, E{B-V)host = 0.37, derived from 
the Lira relation is at odds with the estimate £'(B-y)ho,st = 0.15 
obtained from the shallow TBC relation. They adopt the lat- 
ter value which, after correction, produces B-V ~ +0.2 for 
SN 2009dc near peak. This is broadly similar to the intrinsic 
color we expect for SN 2007if. We note that Yamanaka et a!} 
(12009) do not appear to correct for the efficiency factor a in 
their derivation of the ^''Ni mass of SN 2009dc; however, they 
also assume a fairly typical 20-day rise time for the SN, when 
the true value (not strongly constrained by pre-maximum ob- 
servations) is probably higher, given the SN's slow decline. 
These two errors have opposite and comparable effects on the 



derived ^''Ni mass, so it turns out that their estimate should in 
fact be accurate, and hence comparable to SN 2003fg. Taken 
together with the low, slowly-evolving 50 km s"' d"' to 
day +17) Si II velocity, these observations are broadly com- 
patible with a shell-structure density profile for the ejecta of 
SN 2009dc. The strong detection of C II A6580 in the neai-- 
maximum-light spectra of SN 2009dc, at comparable veloc- 
ity to Si II, suggests either that a larger fraction of carbon 
remained unburned than in SN 2007if, or that a pre-existing 
carbon/oxygen envelope may have been mixed into the ejecta 
after burning was complete. 

The fact that the Lira relation appears to overcorrect the red- 
dening for SN 2007if and SN 2009dc makes the interpretation 
of SN 2006gz problematic. The supposed high luminosity and 
mass of SN 2006gz rests on the Lira relation, corroborated by 
a low-significance detection of Na I and the assumption of a 
steep slope for the TBC relation, and the assumption Rv = 3.l. 
The shallow TBC relation (adopted for SN 2009dc) gives 
E(B -V)host = 0.05, which would make SN 2006gz somewhat 
brighter and redder than average, but not unreasonably so for 
normal SNe la. 

To further examine this question we use the measurements 
compiled in Table^to more closely examine the properties of 
this population under different extinction treatments. To aid 
in this comparison we have remeasured ^^(Na 1 D) from the 
host of SN 2006gz using the same approach as for SN 2007if, 
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Fig. 13. — Constraints on progenitor model masses for SN 2007if, assuming a ■ 



possible central densities. Probability density is normalized to that for the fiducial central density of pc = 3 x lO' g cm ^ 



68% (red), 90%, 95%, 99%, and 99.7% CL (blue). Top left: p,. = 1 x lO" g cm"^ . Top right: pc 
Bottom right: p^ = 3 x lO' g cm"^ 

resulting in a significant (5 x ) reduction in the uncertainty. By 
examining My, Bmax- Vmax, and the color excess, E{B-V)urd, 
with respect to the Lira relation we find that the properties 
of this population are most uniform when the shallow TBC 
relation is used to correct for host-galaxy extinction. The re- 
sulting RMS in Mv is only 0.40 mag, the RMS in B„,ax-Vma.x 
is only 0.07 mag, and the Lira color excess has a remark- 
ably small RMS of only 0.03 mag with this treatment. Use 
of the steep TBC relation doubles or triples the RMS for each 
of these quantities, suggesting that it is disfavored by these 
data. Using the shallow TBC ^^(Na I D) extinction correc- 
tion places SN 2006gz at My = -19.38, well below the popula- 
tion mean of {My) = -20.04. (In fact, without SN 2006gz the 
mean brightness of the population increases to (My) = -20.26 
and the RMS is only 0.16 mag.) 

Other indirect evidence also points to a lower mass and in- 
trinsic luminosity for SN 2006gz. Although the decay time 
of the SN 2006gz light curve is long, the rise time is short, 
lowering the ^^Ni mass estimate relative to what one might 
expect given the event's Am]s( B). Late-time spec troscopy 
and photometry of SN 2006gz dMaeda et al.\ 120091) do not 
provide strong evidence for a large mass of ^^Ni in the ex- 
plosion. The Si II velocity at maximum is typical of nor- 
mal SNe la, which suggests neither a high binding energy 
nor a very massive shell. While the plateau in the Si II ve- 
locity for SN 2006gz supports the existence of an envelope. 



1.3 (see e.g. upper left panel of Figure [TOt for different 

Colored regions are 



: 3 X 10** g cm"\ Bottom left: p^ = 1 x 10** g cm" 



as iHicken et al.\ (120071) suggest, the high velocity of the sup- 
posed plateau (13000 km s~') and its short duration, as well 
as the blue color of the SN relative to the SN 2003fg-Uke 
SNe, argue for an envelope with mass only about 3% of the 
white dwarf mass (0.04 Mq for Mwd = Mch)- Thus, while 
SN 2006gz might have been a DD merger and does show 
some evidence for a low-mass shell, it seems to be differ- 
ent from the red, unambiguously overluminous SN 2003fg- 
like SNe. If the short ri se times of the PDD models of 
iHoflich & Khohkiovl (11996) relative to other explosion mech- 
anisms are realized in nature, SN 2006gz may well have been 
a Chandrasekhar-mass PDD; the models PDD9 and PDD la 
compare favorably. 

In short, comparison of SN 2007if with other events sug- 
gests the following about SN 2003fg-like SNe (excluding 
SN 2006gz): 

1. The tamped detonation scenario, in which the high- 
velocity ejecta are decelerated by an envelope, provides 
a cleaner explanation for the low photospheric veloci- 
ties than a high binding energy alone, and it does so in 
the context of a nearly spherical explosion. 

2. The presence of a shell also explains the long rise times 
and red colors of these events. 

3. While the late-time color evolution of SN 2003fg-like 
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TABLE 4 

Comparison of super-Chandrasekhar-mass SN Ia candidates 



Attribute 


SN 2003fg 


SN 2006gz 


SN 2007if 


SN 2009dc 


Population 












Statistics 


mv 


20.43 ±0.05 


15.99±0.0l'' 


17.26 ±0.03 


15.20±0.16= 




Galactic E(B-V) 


0.013 


0.063 


0.079 


0.071 




Distance Modulus 


40.430 ± 0.005^ 


35.025 ±0.060'' 


37.569 ±0.028 


34.859 ±0.051'' 




Mv 


_2n 04 + 05^ 


_1Q 23 -L n offi 


_2n 3Q + 04 


-19 88 + 17'^ 






+0.14 


+0.03'' 


+0.15 


+0.21 




vsi II (kms"') 


~ 8000 


~ 13000 


~9000 


~ 8000 


9500 (2100) 


Am,5(B) 


0.94"= 


0.69 ±0.04'' 


0.71 ±0.06 


0.65 ±0.03'^ 


0.75 (0.11) 


C II 6580 persistance 


not detected 


+10 (< +11) days 


+7 « +10) days 


+6« +18) days 


+8 (2) days 


Host features 


emission lines; 


Scd spiral; no local 


emission lines; 


SO; no star 






in tidal feature 


emission excess'' 


dwarf 


formation?'' 




Host EW(Na I D) 


< 1.7 A 95% CL 


0.33 ±0.03 A 


<0.12 A95%CL 


1.0 A" 




E{B-V)ur. 


0.21 ±0.04 


0.18 ±0.05 


0.18 ±0.04 


0.37 ±0.08 





Quantities Connected For Host Extinction Using Na I D 

Shallow Steep Shallow Steep Shallow Steep Shallow Steep Shallow Steep 

Con-ectedMv -20.04 -20.04 -19.38 -19.74 -20.39 -20.39 -20.35 -21.43 -20.04(0.40) -20.40(0.64) 

Corrected B„„-V,„o, +0.14 +0.14 -0.02 -0.14 +0.15 +0.15 +0.06 -0.29 +0.08(0.07) -0.03 (0.19) 

Corrected £(B-V)L,ra +0.20 +0.20 +0.13 +0.02 +0.18 +0.18 +0.22 -0.13 +0.18 (0.03) +0.07 (0.13) 



"'Howell et aKHMS) 

^ Hicken et al. 1 2007), with connected distance modulus uncertainty. 
" lYamanaka et aU aOO% 

^ There is no H« emission remaining in spectra of SN 2006gz after subtraction of an interpolated background. 
" Inspection of SDSS images shows uniformly red color. 

' For SN 2003fg and SN 2007if the most probable value, EW(Na I D) = 0, was used 



SNe appears to have a slope similar to the Lira relation 
for normal SNe la, the assumption that the intrinsic col- 
ors of these events at late times are similar to normal 
SNe la probably leads to an overestimate of the redden- 
ing. 

More detailed analysis for these three SNe can be used to de- 
termine likely values for the progenitor masses, allowing for 
the possibility of shells, and (for tamped detonations) to con- 
strain the mass of the envelope. Such an exercise on a repre- 
sentative sample of SNe may be helpful in refining theoretical 
models of the double-degenerate merger process. 

Finally, we note that as an extremely faint system, the host 
of SN 2007if is most likely of low-mass and low-metallicity. 
The tidal tail in which SN 2003fg occurred may have simi- 
lar properties. In contrast, SN 2006gz occurred in a luminous 
spiral and SN 2009dc occurred in an apparently quiesent ellip- 
tical galaxy. Thus, it remains unclear what role the progenitor 
metallicity may play in producing such events. 

7. CONCLUSIONS 

Our observations of SN 2007if provide detailed observa- 
tions of the evolution of a candidate super-Chandrasekhar- 
mass SN la event, allowing new constraints on their progeni- 
tors. SN 2007if is the brightest SN la yet discovered, with an 
inferred ""^Ni mass of 1.6 ±0.1 Mq; the near-lR light curves 
demonstrate evidence for a large fraction of iron-peak ele- 
ments in the ejecta, some of which may have been mixed 
into the shell. The slow Si II velocity evolution near maxi- 
mum is inconsistent with "normal" SN la photospheric evo- 
lution of expanding ejecta with a nearly exponential density 
structure, and is more readily interpreted as evidence for an 



overdense shell in the ejecta. Further evidence for an un- 
usual density structure in the ejecta comes from the late-time 
bolometric light curve, which implies a higher gamma-ray es- 
cape fraction than one might expect for a progenitor with nor- 
mal photospheric evolution. Interpreting the observations in 
the context of a tamped detonation model representative of 
a double-degenerate merg er (as in the DET2ENVN models 
of iKhokhlov. Miiller & H5 flich 20063, we use the SN 2007if 
bolometric light curve to establish the first constraint on the 
total mass for a super-Chandrasekhar-mass SN la progeni- 
tor. Our mass estimate should strictly be construed as a lower 
bound, since reddening by dust and a low central density will 
both result in a higher mass; however, extreme compositions 
and unrealistically high masses, at the edge of the allowed 
parameter space, result if either the reddening or the central 
density are ver y dif ferent from those used in our fiducial anal- 
ysis in Section |53] 

Better models of the progenitors and explosions of 
SN 2003fg-like SNe la are urgently needed, mainly because 
of the theoretical l imits on the existence of supermassive 
white dwarfs ( Pirol l2008t IChen & Lil |2009|) which provide 
possible initial conditions for a SN 2003fg-like explosion. In 
our work, we nevertheless assume the existence of a central 
su permassive white dwarf and use the binding energy formula 
of lYoon & Langed ( l2005h . This formula has been validated 
only for M < 2.0 Mq, but other works have performed sim- 
ilar extrapolations (e.g.. lJefferv. Branch. & Ba ron 2006) sim- 
ply because no appropriate alternative exists for modeling hy- 
pothetical explosions of supermassive white dwarfs. 

Further uncertainty relates to the properties of the 
carbon/oxygen envelope. The envelope considered in 
DET2ENVN and in this work is relatively small in extent 
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(R ^ 10 cm), so the shock interaction which produces the 
shell can be approximated as occurring instantaneously, with 
the ejecta in homologous expansion thereafter Our simple 
modeling suggests that a large fraction of the kinetic energy 
(up to 20% of the initial value, or a few 10^" erg) must be 
dissipated in the transition from explosion to the final state 
with fully developed shell. If the envelope is more diffuse, 
the shock interaction and conversion of kinetic energy would 
proceed over a longer period of time, with some of the energy 
released at optical wavelengths. The interaction could pro- 
duce a "pseudo-continuum" similar to those seen in S Ne 11a 
jHamuvefa/. 2003; Wood-Vasey ef a/. 2004; Alderin g et al.\ 
l2006L.Prieto et aLlQQl) . without the telltale H and He emis- 
sion, but possibly including C and O emission, toplight- 
ing absorption from those elements in the SN atmosphere 
(jBranch et al. 2000). This could explain the blue, relatively 
featureless continuum seen in SN 2007if at -9 days. It would 
also contribute to the long rise time and high luminosity of 
the event, thus inflating our derived ^^Ni mass. Although we 
find this scenario unlikely, very early-phase observations of 
future events similar to SN 2007if, along with theoretical cal- 
culations of synthetic spectra of such an interaction, will be 
needed in order to constrain it. However, even if our ^^Ni mass 
is too high, the uniformity of the optical spectra at late times 
provide some evidence that the interaction produces negligi- 
ble late-time luminosity. The lower limit on the total ejected 
mass from the late-time bolometric light curve should there- 
fore remain secure, since the anticorrelation shown in Table |3] 
and Figure [TT] requires higher total masses for models with 
lower ^^Ni mass. 

More definitive statements about the progenitor mass and 
observational properties await detailed radiative transfer cal- 
culations. These will require more accurate density structures 
of super-Chandrasekhar-mass progenitors, starting from ap- 
propriate explosion models — which may in turn require ini- 
tial conditions appropriately chosen, e.g., from simulations of 
the dynamics of double-degenerate inspiral events. The theo- 
retical uncertainty on the parameter a relating the radioactiv- 
ity and bolometric luminosities at maximum light presented 
another limiting factor, in light of the degeneracy between 
the ^^Ni mass and the total mass for a particular bolometric 
light curve. Clearer predictions for the rise time and intrinsic 
color of super-Chandrasekhar-mass SNe la, with and without 
envelopes or shell structure, would also be helpful in under- 
standing these events. 

A larger sample of SN 2003fg-like SNe is also needed to 
assess whether velocity plateaus truly are generic in these ob- 
jects; this will require fairly intensive spectroscopic follow-up 
as quickly after explosion as possible. Our extensive spec- 
tral time series, including one of the earliest (the earliest if 
SN 2006gz is not included in this subclass) known spectrum 
of a SN 2003fg-like SN la (-9 days), should aid in the prompt 
identification of these events for future follow-up, especially 
in cases where redshifts of faint host galaxies are difficult to 
obtain. Optical/NIR light curves out to 100 days past maxi- 
mum light would be useful to confirm the high mass in these 
explosions. NIR observations could be particularly useful for 
reducing the uncertainty in the dust extinction (and thereby 
testing the Lira relation) for SN 2003fg-like SNe, and for 
determining the distribution of iron-pea k elements for these 
objects. Modeling of the photospheric dMazzali e f a/."2007. 
2008; Gal- Yam er aZ. 20 09) and nebular (e.g.. iMazzali e t al. 
1997l;lMaeda et a/.H2009l) spectral time series would be useful 



to provide independent constraints on the Ni mass, the to- 
tal ejected mass, and the specific products of nuclear burning 
in super-Chandrasekhar-mass SNe, allowing us to learn more 
about the details of the explosion. 

These detailed results will allow us to begin to connect 
SN 2003fg-like SNe to potential counterparts among "nor- 
mal" SNe la. Construction of late-time bolometric light 
curves for a represe ntative sample of le ss luminous velocity- 
plateau and/or LVG dBenetti et a/.ll2005l) SNe la could enable 
new measurements of the ejected masses and shell mass frac- 
tions for these objects, presenting strong new constraints on 
their progenitors. Comparison of these SNe with SN 2003fg- 
like SNe could also enable observational methods to isolate or 
"tag" a subclass demandin g different treatment when measur- 
ing luminosity distances. Oui mbv. Hoflich & Wheeled ( l2007h 
caution that shell-structure SNe should be treated differently 
from normal SNe la in cosmology fits, having intrinsically 
redder c olors and vio l ating the typical width-luminosity re- 
lations. iHowell et al\ (120061) give a similar warning relating 
to a possible class of super-Chandrasekhar-mass progenitors 
with smaller inferred masses than SN 2003fg. Limiting the 
impact of such scenarios is crucial for the success of future 
(Stage 111 & IV) dark energy experiments. Fortunately, with 
the addition of SN 2007if, the properties of this sub-class are 
becoming more clear, and their clear spectral peculiarity es- 
tablishes a means to recognize such events. 
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